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Abstract 

We present a multiscale asymptotic framework for the analysis of the macroscopic behaviour 
of periodic two-material composites with high contrast in a finite-strain setting. Our derivation 
starts with the geometrically nonlinear description of a composite consisting of a stiff material 
matrix and soft, periodically distributed inclusions, where the contrast between the components 
is coupled with the period e. We assume that the deformation field is a rigid-body motion 
on some part of the boundary of the composite, and consider various loading regimes, which 
are controlled by an additional parameter A^. We show that for particular regimes the elastic 
displacements are small throughout the material as e — >■ and that the deformation gradients in 
the stiff material are situated in the vicinity of the rigid-body motion specified by the boundary 
condition. This allows to replace, in the homogenisation limit, the nonlinear material law of the 
stiff component by its linearised version. The behaviour of the deformation gradients on the soft 
component depends on the asymptotics of A,; . In particular, in the case A^ ~ 1 they are shown 
to be of order one, and the related contribution to the homogenised energy is expressed via the 
quasiconvex envelope of the original stored-energy density defined on suitably rescaled gradient 
fields. As a main result we derive, in the spirit of F-convergence, a limit energy functional and 
establish a precise two-scale expansion for minimising sequences. 

1 Introduction 

In the present work we continue the multiscale theme initiated in |14j and study the homogenisation 
of two-material composites with high contrast. We study geometrically nonlinear elastic composites 
with periodically distributed soft inclusions. For applied loads that satisfy a certain smallness 
condition we derive in the spirit of T-convergence an effective two-scale model that — depending on 
the scaling of the load — is either a homogenised, linearised functional, or a homogenised, partially 
linearised energy functional that still allows for large strains in the soft inclusions. In the former case 
we prove that minimiscrs admit a specific two-scale expansion that is determined by the limiting 
model. 

While our work is one of the few papers, along with |14| . [6], that treat the fully nonlinear 
high-contrast case, there has been a significant body of literature during the last decade devoted to 
the mathematical analysis of phenomena associated with, or modelled by, a high degree of contrast 
between the properties of the materials constituting a composite, in the linearised setting. The first 
contributions in this direction are due to Zhikov |39j . and Bouchitte and Felbacq [S], following an 
earlier paper by Allaire [1] and the collection of papers by Hornung et al. [24] (see also the references 
therein), where the special role of high-contrast elliptic PDE was pointed out albeit not studied in 
detail. These works demonstrated that the behaviour of the field variable in such models is of a 
two-scale type "in the very limit" as the period of the composite goes to zero. They also noticed 
that the spectrum of such materials has a band-gap structure and indicated how this fact could be 
exploited for high-resolution imaging and cloaking. It has since been an adopted approach to the 
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theoretical construction of "negative refraction" media, or more generally "metamaterials" , which 
is now a hugely popular area of research in physics (see e.g |34j and references therein). On the 
analytical side, a number of further works followed, in particular [3], [4], [9], [10], [11], [13], [26] . 
|12| . |36| . where various consequences of high contrast (or, mathematically speaking, the property 
of non-uniform cllipticity) in the underlying equations have been explored. Among these are the 
"non-locality" and "micro-torsion" effects in materials with high-contrast inclusions in the shape 
of fibres extending in one or more directions, the "partial band-gap" wave propagation due to the 
high degree of anisotropy of one of the constituent media, and the localisation of energy in high- 
contrast media with a defect ("photonic crystal fibres"), all of which can be thought of as examples 
of "non-standard" , or "non-classical" , behaviour in composites, which is not available in the usual 
moderate-contrast materials. Our goal is to develop a rigorous high contrast theory in a context 
where the underlying nonlinearity is an essential feature of the material description, such as that of 
finite elasticity. In particular, here we focus on the case when the deformation gradients induced by 
the applied forces are situated near the group of rotations SO (3) in the main "stiff" material and 
arc bounded (while being possibly a finite distance away from S0(3)) in the "soft" material of the 
inclusions. This kind of nonlinearity necessitates the development of two-scale homogenisation in 
combination with a rigidity statement in the spirit of Friesecke, James and Miiller [23], which we 
carry out in the present work. 

We next give a more detailed outline of the present article. The starting point of our analysis is 
the elastic energy functional 



where C R'' (for dimension d > 2) is a reference domain occupied by the undcformed elastic 
body, M : — J- R'' is the deformation map, and is a (non-convex) "storcd-energy" density 
that describes the elastic properties of a two-material composite. In the present work, e is a small 
positive parameter that represents the length-scale of the material heterogeneity {i.e. the period 
of the composite) as well as the contrast in the elastic properties of the two material components. 
We consider geometrically non-linear materials and assume that both materials are stress-free in the 
reference configuration. More precisely, we consider stored-energy densities that are frame- indifferent 
and have a single energy well at SO{d) with quadratic growth (see Definition [2T] for the precise list 
of assumptions). 

We consider periodic composites consisting of a "stiff" material matrix (occupying a connected 
subdomain 57^ C il) and e-pcriodically distributed "soft" inclusions (occupying the set = i7 \ 17^ , 
which is assumed to be disconnected). We suppose that on a given part T C dil of the specimen 
boundary the elastic body is clamped, i.e. u(x) = x, and that in its interior the elastic body is 
subjected to a body force. We are interested in the situation where the force density is small relative 
to the strength of the stiff material, so that on the stiff component ftl the deformation is close to 
the identity map. However, due to the high contrast between the two materials, large strains may 
still occur on the soft component. The goal of this paper is to derive rigorously an effective model 
describing the overall behaviour of such a composite on the macroscopic scale. Our approach is 
based on the idea of "splitting" the behaviour of the displacement gradients on the soft and stiff 
components, deriving a priori estimates for these gradients, and rescaling them appropriately. 

To make this precise, suppose that the stored-energy density of the composite is of the form 



where Xe a-nd Wq (resp. xl ^^'^ Wi) are the e-periodic set indicator function and the stored-energy 
density of the "soft" (resp. "stiff") component (so that Xe + Xe ~ l)j ^^d 7 > is a parameter that 
measures the contrast between the materials occupying the two components. Thus, represents a 
high-contrast version of the "classical" stored-energy function {x,F) 1— >■ Wo(F)x'^{x) + Wi{F)xl{x). 
We study the asymptotic behaviour as e of the energy 




W,{x,F) -.^ Wo{F)xUx)+e-^''Wi{F)xl{x), 




(1.1) 
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subject to the boundary condition u{x) = x on F. As a main assumption, we suppose that the loads 
(Is) C L^(il)'' \ {0} satisfy the following "smallness" condition: 

limsupAe<oo, where K ■= £''\\ls\\L^{n) + £\\le\\L'^{no)- (1-2) 

As part of our main result wc prove that elements of minimising sequences (ug) for admit a 
representation in the form 

u,{x) =x + Xe{eg°{x) + e'^gli^)) with {glgl) e H^{^}°)xH^^„{n), 

where Hp ^{il) := {u e H^{n) : u\r — 0} and iJQ(f2") denotes the closure in H^{il) of smooth 
functions with support in fi^. The pair [g^.g]) is interpreted as a "splitting" of the displace- 
ment Ue — id where the soft-part contribution g^ and stiff-part contribution g^ are scaled differently, 
according to the scaling of the contrast. The heart of our analysis is an a priori estimate (see The- 
orem |3]4]). It is based on a geometric rigidity estimate for perforated domains (see Proposition 15. 2p . 
which is a specially adapted version of the geometric rigidity estimate derived in the work |23] by 
Friesecke, James and Miiller. Our a priori estimate shows that the sequences (g°), (eVg°), (gl) and 
(Vgl) are bounded in i^(il). As a consequence, we have, up to a subsequence, the weak two-scale 
convergence 

5° ^ y), eVg° - V,g°(a;, y), \7gl - \7g\x) + \7y^{x, y), 

where g" G L^{n; H^{Yo)), g^ G H^.oi^) and G L^{n;H^{Y)). Here Y := [0,1)^ the set Yq d Y 
represents the soft component, and H^{Y) is the Sobolev space of periodic iJ^-functions. The main 
result of the present paper is that the limit ((/"jt?^) is naturally related to a functional of the form 

A(ff°,5')-/ l-{g"+g')dxdy, (1.3) 

where A is a "homogenised" energy, which wc rigorously derive from A^ in the spirit of two-scale 
F-convergence. We can extract certain information on the precise form of the expansion ([T]) from 
the minimisers of (|1.3p . The quality of this information and the functional A itself depend on the 
scaling properties of the applied load Z^. In particular, we distinguish the cases 

Ae ^ 1 as e — > 0, and — > as e 0. 

We first discuss the regime A^ — > 0. In this case the behaviour of both components is shown to be 
accurately captured by the linearisation of the storcd-energy densities: we have (Theorem 14. 3p 

A(5«,gi) = A-°^"(g0,5i) := (^^^ (v,.g°(2:, y)) dy + QY'^^ {w g\x))^ dx. 

Here Qi denotes the quadratic form from the expansion of Wi at the identity, i = 0, 1; further, 
Qhom result of "perforated" homogenisation of Qi. By appealing to the convexity of the limit 

problem, we show that the splittings {g^,g\) associated to minimising sequences are strongly two- 
scale convergent to the minimiser of the limit problem, say {g'^,gl). Broadly speaking, this means 
that any minimising sequence can be approximated by the expansion 

ue X + X,(eg°{x,x/e) + e'-'iglix) + eiP^)'^ (1.4) 

with a corrector -0* that is determined by the cell problem associated to Q\°™ (Proposition 14. 4[) . 

In the regime Ag ~ 1 only the material behaviour on the stiff component is linearised and A takes 
the form 

A{g^,g') ^ Ai-8''(g°,.gi) := (^^ Wo[l + Vyg%x,y)) dy + g'r"(Vgi(a;))) dx. 
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In this case our result is weaker: we prove that there exists a minimising sequence for £^ such 
that 



Due to the non-convcxity of the storcd-cncrgy density, we expect that a general minimising sequence 
may oscillate on scales larger than e and therefore cannot be captured by the two-scale ansatz above. 
However, we prove the convergence of the sequence of infima of to the minimum of the limit 
problem (Theorem 14. 7p . 

In the rigorous formulation and the proof of the above results we use a number of mathematical 
tools that have been developed in the last two decades in the areas of calculus of variations and 
homogenisation. Among these are the two-scale convergence introduced in [35] , [I] ) and the periodic 
unfolding (see |17] and references therein), the basic information about which is given in Appendix 
A. These provide an appropriate vehicle for the description of the two-scale structure of low-energy 
sequences in our setup, which at the heuristic level is represented by the expansions (|1.4p and (|1.5p . 
Further, in order to formulate a suitable two-scale compactness principle for such sequences (see 
Proposition 13.101 below), we make use of an extension theorem for periodic perforated domains 
(see [33], [7]) in combination with a geometric rigidity statement (see [23]). Finally, we exploit 
substantially the mathematical framework of F-convergence (see [TH] and references therein), which 
allows us to formulate a rigorous statement about the convergence of (rescaled) energies, see Theorem 
4.1 and Theorem 4.5. 

The paper is organised as follows. In Section [5] wc list our assumptions on the material law and 
the geometry of the composite microstructure. In Section [3] we introduce a scaled splitting of the 
displacement that is adapted to the composite geometry and to the contrast in the clastic stiffnesses 
of the two components. Further, we establish an a priori estimate for the splitting. In Section 13.31 
we prove weak two-scale compactness for low-energy sequences, identify the corresponding limits 
and present an approximation scheme for the limits. Section [4] is concerned with the F-convergence 
analysis, where the regimes ~ 1 and arc considered separately. Finally, in Section [5] we 

present our adaptation of the geometric rigidity estimate to perforated domains. Some relevant facts 
about two-scale convergence and some auxiliary results are reviewed in the Appendix. 

We use the following notation (see Section [2] for the details): 

• d > 2 is the dimension; 

• C M'' is an open, bounded, connected Lipschitz domain; 

• e is a small positive parameter equal to the period of the composite; 

• (resp. ill) ^oft (resp. stiff) component; Xe (resp. xl) is the associated indicator 



• Y := [0, 1)'' denotes the reference cell of periodicity; Yq CC Y represents the (re-scaled) soft 
component; Yi := Y\Yo represents the (re-scaled) stiff component; is the indicator function 

for Yi, i = 0, 1; 

• LP{n;R'^), H^{n;R'^), p(17;R'') and H^{Yo;R'^) denote the usual function spaces of vector 
valued, p-summable functions, square-summable Sobolev functions, the closure in _ff^(f2;R'') 
of functions from C^(fi;R'') vanishing in a neighbourhood of F C dfl and the closure of 
Cq{Yo;R'^) ill H^{Yo;R'^), respectively; H^(Y;R'^) denotes the space of F-periodic maps in 
Hl-^^(R'^;R'^). The function spaces Geifl) and GrfliflxY) arc defined in Definition 13.11 and in 
(|3.18|) below. If it is clear from the context we drop the target space and sometimes even the 
domain of definition in the notation, e.g. we write L^(fJ) or even for L'^{il;R'^) as well as 
for L2(f7;Krfxrf). 

2 

• We write ^ and — > for the weak and strong convergence in a Banach space. We write ^ and 

for the weak and strong two-scale convergence (see Definition IA.3P , and we write and 

sp2 

— >■ for the weak and strong two-scale convergence in the "splitting" sense, as introduced in 
Definition im 




(1.5) 



function; 
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• We write f{h) = 0{h°') whenever limsup^_^Q < oo, and f{h) = o{h°') whenever 
Hni,.^o/j-"|/(MI =0; 

• For real ((ixd)-matrices F = (F^) and G = (Gij) we write F : G :^ PijGij for their inner 
product, symF := (F^ + F)/2 for the symmetric part and skwF := F — symF for the skew 
symmetric part of F\ further, we denote by S0(o?) the set of rotations in W^. 

• Unless stated otherwise, we write C for a generic positive constant which is independent of e 
but can change its value even within a single chain of inequalities. 

2 Setting of the problem 

We first introduce our geometric assumptions on the sets occupied by the soft and stiff materials. 
Recall that in a reference configuration the composite occupies a domain 17, which is assumed to be 
an open, bounded, connected set in M'', d > 2, with Lipschitz boundary. In order to describe the 
microstructure of the composite, denote by Y :— [0, 1)'' the unit cell of periodicity and suppose that 
Yq C y is an open set with Lipschitz boundary. We assume that Iq is compactly contained in 
i.e. A\st{YQ^dY) > 0, and set Yi ■.— Y\ Yq. Next, for £ > we decompose f2 into a "stiff" matrix 
Vl\ and "soft" inclusions $7^, as follows: 

^l.= \J{e{i + Y^) : CGZ^£(^ + F)cf7}, nl:=^l\^l. 

By definition ilj is connected and \ 517" = 912 for all e. Furthermore, 12^ is a disconnected 
set consisting of small inclusions of size ^ e. Note that the soft inclusions are eF-periodically 
distributed in the interior of VL and that they have a "safety distance" to the boundary 9f2, i.e. if a 
cell + Y) intersects dVL, then + y) n 12 is occupied only by the stiff material. We denote by 
X\ the characteristic functions of the sets f2^, i = 0, 1. 

For given parameter values e > 0, 7 > and a load £ L^(r2)'', let the integral functional 
£^ : (17; W^) M be given by (jl.ip . As is mentioned in the introduction, the parameter 7 measures 
the contrast between the soft and stiff component, and e is a small parameter that determines the 
period of the composite. In the limit e — > we expect the energy to be (partially) linearised; 
therefore, it is convenient to rewrite the elastic energy part of £e{u) as a function of the displacement 
ip(x) := u{x) — X : 

A,((^) (Wo{I + Vip)x° + e-'Wi(/ + V(p)x^) dx. 

Clearly, with this convention we have 

£e{^<^ + y^) ~ ^eiv) - / le-ifdx, 

where id(a;) := x denotes the identity map in M''. We are interested in the asymptotic behaviour 
of Ag along low-energy sequences of displacements that satisfy the clamped boundary condition on 
some part of the specimen boundary. More precisely, let F be a measurable subset of dfl with a 
positive {n — l)-dimcnsional Hausdorff measure. In what follows we set 

w^:^{n) ■.= {^&W''^{n) : ^Ir-O}. 

To avoid technical difficulties we assume throughout the article that F is regular enough, so that 
Wj!;'^(12) is dense in o(^)- The boundary condition for the deformation u and the displacement 
(fi then reads 

u - id = e i/p o(12). 

Unless stated otherwise, (e) denotes an arbitrary, but fixed vanishing sequence of positive num- 
bers (ej)jgN- Quantities {g^) indexed s are understood as sequences of the form {g^^)j^jq] in the 
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same spirit e denotes a generic element Ej of the sequence (Ej ). Since the convergence results and 
estimates are independent of the specific choice of the vanishing sequence (ej), they are generalised 
in a straightforward way to the setting with a continuous parameter e. 

Next we list our assumptions on the material law. In addition to the usual properties assumed in 
finite elasticity (such as frame- indifference, growth and regularity), we require that both materials 
be stress-free in the reference configuration (i.e. the reference configuration is a "natural state"). 
Unless stated otherwise we assume that both energy densities Wq and Wi belong to the class Wg, 
> 0, as defined below. 

Definition 2.1. The class We, d > 0, consists of all continuous maps W : M'*^'^ — > [0, oo) with the 
following properties: 

• W is frame-indifferent, i.e. W{RF) = W{F) for any R e SO{d) and any F e E''^'*; 

• I is a natural state, i.e. W{I) = minVF = 0; 

• W is non- degenerate, i.e. for any F E R'^^'' 

W{F) > 9dist^{F,S0{d)); 

• W admits a quadratic expansion at the identity, i.e. there exists a non-negative quadratic form 
Q onR'^'"^ such that 

W{I + G) ^QiG)+r{G) (2.6) 

for any G G M''^'*, where 

r{G) = o{\G\^) as |G| ^ 0. 

Remark 2.2. Let W he of class We o,nd let Q be the associated quadratic form. Then Q vanishes 
for skew- symmetric matrices and is positive definite ( with a constant that only depends on d and 9 ) 
on the subspace of symmetric matrices. Thus the integral J^-^ Q (V (p) dx is interpreted as the elastic 
energy of the displacement ip e H^{n) in the framework of linearised elasticity. 

Remark 2.3. For the a priori estimates that we derive in Section\^the expansion condition ()2.6p 
can be dropped. In order to simplify some reasoning for the regime \e ^ 1 we additionally assume 
in Section \4^.S\ that Wq has quadratic growth at infinity and is quasiconvex. 

Remark 2.4. Under the requirement of material frame- indifference, i.e. W{RF) = W{F) for any 
{d X d)-matrix F, such that dctF > 0, and any rotation R, we have W{F) = W{F'^ F) for some 
function W, and the first Piola-Kirchhoff stress tensor T at a point where Vw ^ F is given by 
T{F) = VM^(F) = 2F\'WiF'^F). 

Remark 2.5. Suppose that W : R'^'^'^ [0, oo) is frame-indifferent, non-degenerate, stress-free at 
the identity I and of class G^ in some neighbourhood of I . Then W is of class We and one has 
Q{G) = D'^W{I){G,G)I2 > G|symG|2 for some G > 0. 

3 An a priori estimate and two-scale analysis for low-energy 
sequences 

We consider the minimisation problem for the energy 

As{ip) ~ / l^-ipdx, where 99 g _ffp Q(f2). (3.7) 
Jn 

Throughout this section we assume that (l^) C L'^{ft)'^ \ {0}, and that Wi, i = 0, 1, arc of class We, 
although we need not require Wi to fulfil the "quadratic expansion" condition (|2.6p . The main result 
of this section is an a priori estimate for sequences with low energy. The estimate is based on a specific 
splitting (combined with a certain scaling) of the displacement tp into a contribution determined by 
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the stiff component and the remainder. Motivated by this, we introduce in Section TS-Sl several notions 
of two-scale convergence that are adapted to the spUtting, prove a two-scale compactness result and 
identify the admissible limits. Furthermore, we present an approximation scheme for the admissible 
two-scale limits. Let us remark that the r-convergence analysis in Section 2] strongly relies on the 
two-scale results derived in the present section. 



3.1 A note on the natural scaling of "near- minimising" sequences 

Let us briefly describe the main idea behind the analysis to follow in Sections[3]and|4l The homogeni- 
sation and asymptotic analysis of p.7p is carried out in the framework of F-convergence introduced 
by De Giorgi (see [IHlIin] ). In order to obtain relevant information about the asymptotic behaviour 
of (|3.7|) by this approach, it is mandatory to carefully select a "natural" notion of convergence for 
the underlying deformation fields. As in [20], the right object to study are displacements rescaled 
in a way that is adapted to the energy functional. However, since we treat composites with high 
contrast, the right choice of the scaling is not obvious. As we describe below, our choice is naturally 
motivated by an a priori estimate of the displacements ifg = — id and a splitting of ip^ adapted 
to the geometry of the composite. 

We next give a rough idea of how the a priori estimate mentioned above is obtained. The 
non-degeneracy of the components directly implies 

II dist(VH„ SO(d))||i.(j,o) + e-'^ll dist(Vwe, SO(d))||i.(i,i) < CA,(^e). 

Since 7 > 0, the second term on the left-hand side is of leading order as £ 0. Thus, we deduce that 
the behaviour of Vug on the stiff component is dominant. Furthermore, we find that Vu^ is close 
to SO((i) on the stiff component, as can be quantified by using a version of the geometric rigidity 
estimate of [23] adapted to perforated domains (see Theorem 15. ip . In order to gain control of the 
contribution of Vit^ on the soft component, we introduce the splitting it^ = id + ipg = id + ip'^ + (pi , 
where (p^ vanishes on ft^ and pi is harmonic in fl^. We are going to prove that 

llVy^^lli^foo) +£-'^||V^i||i.(,,) < CA,(^,), 

and that the linear, load-related part of the functional (|3.7|) is estimated as 



L ■ ipf, dx 



n 



< 



CA,(||V¥.^|U.(oo)+£-''||V^i|U.(n)) 



where = e'^||Ze||^2(f2') -I- £||^e||L2(-Q0). In conclusion, we obtain the a priori estimate 

||V¥.°|U2(no) +£-^||V^i|U2(a) < CAe, 
provided the sequence of displacements ip^ satisfies 

limsup —r [A£:((y9e) — / l^-ip^dx] < 00. 
6^0 Xj \ Jn J 

In this paper we consider only the cases when A^ is either of order one or converges to zero as 
£ ^ 0. In view of the previous estimate, it is natural to consider the rescaled sequences A~^</5° and 
X^^e^'^'pl. Indeed, the rescaled contribution on the stiff component, \~^e~^(pl, is of order one and 
entails a nontrivial limit. However for the soft component, by a improved Poincare inequality, see 
Lemma [3^ we have 1 1(^°| 1^2(^0) < C£||V(/J°||l2(qo), and we deduce that Aj-^i/j^ converges to zero 
in the L^-norm, although its contribution to the elastic energy is of the same order as that for the 
stiff component. For this reason we consider the scaling \~^e~^tp^g instead. In conclusion, we are 
going to show the following: as £ — >■ 0, the infima of the energy functional p.7p are of order A^, the 
related minimising sequences can be represented as p^ — Xe{£g^ + s'^ gl), where vanishes on fi^, 
gl is harmonic in 57^, and the sequences 5^, £Vgg, gl, Vg^ are bounded in the L^-norm. Therefore, 
we consider the F- limit of the rescaled energy sequence A~^ (^A^{ip^) — J^le ' 'fie dx) with respect to 
the two-scale convergence of the sequences 5°, gl- Let us point out that any other scaling either 
leads to a loss of compactness or to trivial limits. Hence, (for A^ < 1) the scaling introduced above 
is indeed the natural one. 
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3.2 The a priori estimate and a splitting of the displacement 

Consider a sequence of displacements {(pe) in with cqui-bounded energies, i.e. A^{(pi,) = 

0(1). Then clearly the sequence {(pe) is bounded in H^{il). Indeed, since the stored-energy function 
is non-degenerate (see Definition 12. ip , it satisfies a quadratic bound below. In the high-contrast 
situation considered in the present article, we can expect more, since the strength of the stiff com- 
ponent increases as £ — > 0. Thus, it is natural to expect that we gain additional control on the 
displacement gradient restricted to the stiff component fi^. To make this precise by means of a 
quantitative estimate, we introduce a "splitting" of the displacement that reflects the difference in 
the behaviour of the soft and stiff components. Let us remark that this splitting will also play a 
crucial role in the asymptotic analysis of the energy as e 0. 

Definition 3.1 (Splitting). For s > set ^.(fJ) := H^{n°)xH^^Q{n). For each ip e i?r_o(^) ^'^sre 
exists a unique pair (17^,(7^) € Ge{^) such that 

if ~ egl + e'' gl in fi, g^ = Q in fJ^, and 

l|V5,'|U2(o) < \\V{gl +0\\mn) for all ( G Hl,{n°). (3.8) 

We refer to {g^,gl) as the e-splitting of p and write p {g^,gl). 

The splitting is well defined due to the following classical extension result for perforated domains: 

Lemma 3.2. (see e.g. \33l Theorem 4-'^])- Let e > 0, g ^ H^{fll). There exists exactly one map 
g G H^{fl) with g = g on fll and 

[ V5 : VCdx = for all ( G i?o(r2°). (3.9) 

Furthermore, there exists a constant C independent of e and g such that 

\\g\\mn)<C\\g\\L^^ni). (3.10) 
l|V5lU=(o) <C||V<7|U2(oi), (3.11) 
||symV5||i,2(n) < C|| sym V5||i2(oi). (3.12) 

Remark 3.3. Although property p.9p is not stated in Theorem 4^.2 of \33f . it is implied by the proof 
of Lemma 4-1 of \33^ . which is utilised by the latter. 

By definition gl is the unique harmonic extension of e~'^ipe\ni to fl, i.e. the map in H^{fl) that 
satisfies gl = s~^ipe in fij and the elliptic equation p.9p . which is the Euler-Lagrange equation asso- 
ciated to the minimisation problem p.8p . The important property here is that the condition allows 
to control the L^-norm of the gradient of the extension by e~^'^'p restricted to the stiff component. 
Since p.9p admits a unique solution, the splitting defined in this way is unique. Furthermore, note 
that the mapping ^{il) B p 1-^ {g°,gl) G Ge{^) is a linear bijection. 

The above splitting construction is tailor-made for the following estimate. 

Theorem 3.4 {A priori estimate). Let ip G o(^) '"^'^ id^Tdl) V- Then 

\\9%W) + W^&^m + WaWW-m ^ ^Ae(^) (3.13) 

for a positive constant C that only depends on the geometry offl,T, Yq and the parameter 9. 

We present the proof of the above theorem in Section 13.41 In a nutshell, the argument is as 
follows: suppose that {(pe) is a sequence with Ae((^e) = 0(1). By virtue of the non-degeneracy of 
the stored-energy functions and due to the high contrast between the soft and stiff component, the 
scaled displacement gradient e~'''V(pe restricted to the stiff component fll is bounded in L^. Then 
by Lemma [3.21 we conclude that ||5£||i/i(f2) = 0(1). Furthermore, by construction the difference 
ipi, — £^ gl (which coincides with eg^) is only supported by the soft component and can be identified 
with a map in i/o(r2°). Since £7^ is a union of small inclusions of size ~ e, it follows (by using a 
suitable version of the Poincare inequality) that ||5el|L2(f2) = 0(1). 
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Corollary 3.5. Suppose that (Ig) satisfies (|1.2p and {(pe) is a sequence in Q{fl) such that 



limsup ( Ae(i^£) - / le-(pedx 



< oo. 



The 



linisupAg A^{ipi.) < oo. 

e->0 



(3.14) 



(3.15) 



Proof. Let {g^,gl) Lp^ and recall that :~ £''\\le\\L'^(iYj + ^||4||L2(of)- From the definition of 
e-splitting we have 

[ k-p,dx^ f{eV,)-gldx+ f (d,)-5e"dx. 
Jo Jn Jn° 

Using the Cauchy-Schwarz inequality, the a priori estimate p.l3p and the definition of yield 

(3.16) 



le ■ ipe dx 



ai < \-' A,(^,) - / l,-p,Ax]+ \- 



For Oe := A^ ^ y'A^^Jp^ we now have 



• ipi; dx 



< C'il + a,), 



and p.isp follows. 

The bound p.l4p is a low-energy condition. It is satisfied, in particular, if 



AeiVe 



le ■ (fie dx < 0, 



□ 



(3.17) 



when we even get the bound Xj'^A^{(pi;) < C for any e, where C is a positive constant which does 
not depend on the choice of a particular sequence (ip^) satisfying p.l7p and only depends on the 
geometry of r2,r, Yq and the parameter 9. Furthermore, the set of sequences of displacements 
ipe G -ffp g(i7) satisfying p.l4p is non-empty: for example, the trivial displacement ipi; = fulfils the 
inequality p.l7p . 

As mentioned above, in the proof of Theorem 13.41 we need the following improved Poincare 
inequality. 

Lemma 3.6. Let the sets fl, 17" and ftl be as in Section\^ For all u G iJp g(fi) we have 

< C{\\Vu\\l2(q_i) + e\\Vu\\L2(^^o)). 
where C denotes a positive constant that only depends on the geometry ofQ,T and Yq. 



Proof. Let u denote the harmonic extension of itj^^i to fl as introduced in Lemma l3.2l Set ^ := u — u 
and := { ^ G eZ'^ : (, + eYq C }. By construction, C is supported in fi^ = U^ez, + ^^o)- For 
C G and y G Yq define V^{y) := C(C + £J/)- We note that G H^{Yo) and Vyl4(y) = £VC(C + ey). 
Therefore, re-scaling and Poincare's inequality yield 

[ \C\' dx ^ e'' f m'dy<Ce''Y. f l^yVil' dy ^ Cs' [ |VCp dx. 



By the identity u = u + and Poincare's inequality applied to u, we get 

ll"llL2(n) < \\u\\LHn) + \\C\\L^{ni) < C'(||ViI||i2(o) + £||VC||L2(fjo)) 

= C(||VM||i2(o) -|-£(||VM||L2(no) + ||Vw||L2(nO))). 

Hence, the result follows from the estimate p. lip . 



□ 
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Remark 3.7. In high- contrast homogenisation problems one often has to work with sequences of 
functions whose gradients are allowed to be large on some periodic components of the domain. In 
this situation the standard Poincare inequality may not be good enough to provide boundedness of 
Ue, and one has to use its improved versions, such as || Vue||/^2(f^i) + ||eVwe||/,2(-Q0) < C. For further 
details we refer to USf . where a "high-contrast" Poincare inequality is proved under the assumptions 
that has zero mean and is not required to be disconnected. Lemma \3.€\ is a simplified version 
of this high-contrast Poincare inequality. 

3.3 Two-scale limits of low-energy sequences 

Let ((pe) C Hy o(f^) have equi-bounded energy {i.e. it satisfies p.l4p ) and let {g^.gl) denote the as- 
sociated e-sphtting. CoroUarv l3.5l and the a priori estimate in Theorem [33] show that the sequences 
((7°), (eV(7°) and (Vg^) are bounded in and therefore they arc also relatively weakly two-scale 
compact. In this paragraph we identify the corresponding two-scale limits (see Proposition I3.10p 
and present a scheme that allows to approximate the admissible limits (see Proposition l3 . 1 3p . These 
limits belong to the space on which we define the limiting homogenised functional in Section |4l 

We start with the definition of the limit space and certain variants of two-scale convergence 
adapted to the splitting introduced in Definition 13.11 see Definition IA.3I for the usual definition of 
two-scale convergence of functions. Set 

xl^ L^(Sl-Hl{Y)) : II \/ytlj -.X/yCdxdy ^0 for all C e L^{n;H^{Yo))\ . (3.18) 

With slight abuse of notation we write (5°,(7^) € Gr,o (resp. -0 G ^r.o), if (3°,. 9^,0) G Gr.o{^ x y) 
(resp. (0,0,V)eer.o(f^xy)). 

Definition 3.8. Let (ip^) be a sequence in q{^) , {g'^ , g^ E Gr.o{^'X-Y)i o,nd let {g^ , gl) denote 
the e-splitting associated to ip^. 

1. We say that a sequence (ipe) weakly two-scale converges in the splitting sense to {g^,g^), and 
write (fie"^ {g°,g^), if 

e\7g°^Vyg° in L^iflxY), g] ^ g^ weakly m H^{n). 

2. We say that {(p^) weakly two-scale converges in the splitting sense to {g'^ , g^ , ip) , and write 

^-^^ (.9*',ff\'/'), -ifVe {g°,g^) and 

Vgl ^ Vg^ + Vyij in L^{nxY). 

3. We say that (ip^) strongly two-scale converges in the splitting sense to (g*^,.?^), and write 
p, (5°,ffi), tf^e [g^g') and eVgO 4 V,gO m L^QxY). 

4- We say that (ip^) strongly two-scale converges in the splitting sense to {g^,g^,ijj), and write 

fe "-^ {g°,g^,ip), if'Pe {g°,g^) and 

^gl ^ + Vy^ in L'^inxY). 

Remark 3.9. Our requirement of two-scale convergence of the scaled gradients eVg^ in the above 
definition implies the two-scale convergence of g^, which can be seen from the Poincare inequality. 
It leads, however, to an equivalent convergence statement for the functionals A^^A^, due to the 
fact that the sequence {eVg^) is L^ -bounded whenever g^ correspond to a sequence (pe for which 
limsup A~^Ag((pj) < 00. 
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The above definition introduces a hierarchy of convergence types: In particular, the following 
implications hold: 



The weak convergence types (p^ {9'^, 9^) and (p^ {9^,9^,4') are naturally related to the energy 
by means of the following compactness result. 

Proposition 3.10. Let (cp^) C Hp and suppose that limsupg_^Q A^^A£((pe) < 00. Then there 
exists a triple {9^,9^,11;) £ Qr.oi^y^Y) such that, up to a subsequence, one has X~^ip^ {9^,9^,11^). 

Proof. Consider the e-splitting {9^,91) «^ A^^cp^. By Theorem l3.4l the sequences (g^), (.gg), (eV^^) 
and (V^g) are bounded in Thus, by Proposition IA.5I there exist {90,91) G Qvfl and ij} G 

2.2(17; iJ^(r)) such that 

.9° A sVgO A V,/, 

.9^1 ^ 9^ weakly in V9I ^ V 9^ + Vj^V', 

for a subsequence that we do not relabel. It remains to show that the function if) satisfies the integral 
identity 

// \Iy4{x,y):VyC{x,y)dxdy = Q for aU C € i2o'(^o)). (3.19) 

J JOxYo 

By invoking a density argument it is sufficient to take smooth test functions in (|3.19p . i.e. C S 
C^(17; Cq (Yq))- Set Ce{x) := eC,{x,x/e) where C is extended in its second variable by zero to Y 
and the resulting function is extended by F-periodicity to M''. Then clearly, € H}^{Vt^) for small 

enough e and VCe 4 VyC in L'^{VlxY). Thus, 

/ V9I : VCs da; ^ // (Vp^ (x) + Vyi,{x, y)) : VyCix, y) dx dy. 
Jn JJnxY^ ' 

as e — )■ 0. By p.9p the left-hand side of (|3.3p is zero, hence it suffices to notice that 

V.9'(.t) : VyC(a;,y)dxdy = / ^9'{^) - [ VyC{x,y)dydx = 

Jn Jy 



nxY 



in view of the fact that C{x, •) e H^{Yo) C H^{Y) and hence ^yCix, y) dy = 0. □ 

The limit (<?'', <?^) captures the relevant behaviour of {(p^) as e — ^ 0. In particular, the map 
resolves as a two-scale limit oscillations that might emerge in the soft component. We are going to 
show in Section 0] that the limiting behaviour of the energy A^{(pi;) can be described as a function 
of (9^,9^). By combining Corollarv 13.51 and Proposition 13 . 101 we arrive at the following statement 

Corollary 3.11. Let {ipg) C q{Q,) be such that ^3.14^ holds, where Ae, le satisfy the condi- 
tion U.l^) on the applied load. Then, up to a subsequence, X'^^cpe (9^,9^) for some {9^,9^) G 
GrflinxY). 

It might seem that the limit {9'^, 9^) of a sequence {Xj^(p^) depends on the definition of the e- 
splitting {91,9'^) associated to X^^ip^. That is if we allow a certain degree of freedom in the splitting 
construction, namely in the way of extending e^'^ X'^^ipsln'^ to fl, but still requiring p.l0p - p.l2p . 
The next corollary shows that {9^ ,9^) is in fact independent of the specific construction and can be 
recovered from Lp^ directly. 
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Corollary 3.12. Under the conditions of Provosition [3.10\ one has, up to a subsequence: 

e-^X-'^,xl^\yi\g' m L^{n), and X-'V^,^Vyg°. 

Proof. Let {g^,gl) Xj^ip^. For the first convergence we use the fact that, up to a subsequence, 
gl — s- g^ strongly in L'^{il) and xl ^ weakly in L^{il). The second convergence follows directly 
from Proposition 13.101 by recalling that e'' = o(l). □ 

Note that the above corollary fully identifies the map g° as well as g^ : indeed, suppose that g° 
and g° in L'^{il; Hq{Yo)) satisfy Wyg'^ = '^yCj^- Since g^ = 5" on VIxOYq, the Poincare inequality 
then immediately yields g^ = (f . 

Finally, we present an approximation scheme for (5", 5^) G Qt.o(^'><Y) that will be used in 
Section 21 for the construction of recovery sequences in the small- load regime. 

Proposition 3.13. Let (g^^g^^ip) e Qr.oi^'xY), and let (cg) be an arbitrary vanishing sequence of 
positive numbers converging to zero. Then there exists a sequence (ipe) *^ ^(fl) such that 

A^Vs (g°,.9\V'), limsupCe||V¥.,|Uoo(n) = 0. (3.20) 

The proposition shows that any pair {g^,g^) G Qr^i^xY) can be approximated by a sequence 
of displacements ((p^) C q{^) with respect to a convergence notion that is even stronger than 
the one introduced in Definition 13.81 We would like to remark that the second condition in p.20p is 
a technical property that is tailor-made for the linearisation of the stored-energy density, which will 
be at the centre of the argument of Section^ Furthermore, in the proof of Proposition we need 
the following continuity result: 

Lemma 3.14. Let 3/2 < p < 2 and [i]^) be a sequence in H^{il) such that 

\\Ve\\w^.p{ni) ^ 0, (3.21) 

/ V?7e:VCd.T = for all C e H^{n"^). (3.22) 

Then (r]^) strongly converges to zero in W^'P{il). 

Proof of Proposition [3A3[ Recall that ip^ {9'^,9^) means 

mL^inxY), eWg"^AWyg° in L^{nxY), (3.23) 
gl g^ weakly in H^ifl), Vgl 4 V.g^ + V.y?A in L^iflxY), (3.24) 

where (g'^^gl) (p^. The proof of the proposition is divided into three steps. First, we present a 
construction for limits {g'^,g^,'ip) that arc subject to additional regularity assumptions. This yields 
a sequence (pe = eg" + e'^gl, where gl does not necessarily satisfy p.Sp . Next, we extend the 
construction to general limits using a diagonalisation argument. Finally, we argue a posteriori that 
the associated e-splitting satisfies p.23p and p.24p . 

Step 1 (Construction for L°°-data). First, we present a construction for uniformly bounded limits. 
Suppose that 

.9" e C,i(17xro), .9^ e Wr;o°°(l^), and ^ c,Cl{n-Cl{Y)). (3.25) 

Recall that {x} € Y denotes the unique vector such that x = [x] + {.t} for some [x] e (see 
Appendix \K\ . Define 

:=g°(x,{a;/£}), gl{x) g\x) + ei:{x,x / e), p^,{x) := e^,{x) + e''gl{x). (3.26) 
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By construction we have g^,gl,(pg G ffpg(fJ); and if e > is sufficiently small, there holds G 
H^ifl'^,)- Furthermore, ^(3^3^ -dSJl are fulfilled with {g°,gl) replaced by ig^,gl). 

Step 2 (Diagonal sequence construction). Next we approximate general limits by employing a diag- 
onal sequence construction. Let {g'^,g-^,ip) be as in the statement of the proposition. By a density 
argument for each /s S N wc can find maps (g^, g^, ipk) that satisfy p.25p and approximate {g^ , g^ , ip) 
so that ^ 

Il5° -5fellL2(n;ffi(yo)) + ||.g^ ^ glWH^n) + Wi^ - i^kllmn-.H^Y)) < ^■ 

Let (g^^), (gk e) ^^'^ iVk.e) be the approximations associated to {g'^, gj.,ipk) by the construction 
(P:^ . Define ' 

Ck,e := \mgt,e) ~ g"\\L- + imeVgt^,) - V,5°IIl^ 

+ II5L - + We{ygle) - + ^v^)\\l^- + Ce||V(pfc,e||L-(0), 

where denotes the unfolding operator (see Appendix IXj) and = L'^{M.'^xY). (Note that we 
tacitly extend the maps defined on ilxY to M.'^xY by zero.) It follows from Step 1 that the unfolded 
quantities in the definition of Ck^e strongly converge as e — > 0, hence 

limsupcfe^e < \\g''k~g°\\L- + \\^yg°k~^vg°\\L^ + \\gl-g^\\L- + \\{^gl+^yiJk)-{yg^+yymL^- < I- 

Hence, limfe_^oo lim£->o Cfc.e = 0. Now, an argument from [5] (see Lemma IA.6P yields a diagonal 
sequence (fc^) C N such that 00 and Cfc^^e — ?► as e — )■ 0. As a consequence, the sequences 

defined by g^ := g^^ ^, gl := gl^ ^ and <pe := tpk,,e satisfy p.23p - p.24p with (5°, 3]) replaced by 

Step 3 (Convergence of the e-splitting). Let {g^,gl) (pe- Note that in general ((g^,gl) might 
differ from {g^,gl), since gl is not necessarily harmonic in $7^. Thus, we have to argue that the pair 
{g^,gl) satisfies p.23p - p.24p . Setting 77^ := gl — gl, it clearly suffices to prove that 

WVeWv^in) = 0(e), V77e -> strongly in L^{n). 

This can be seen as follows. Since vanishes on the stiff component $7^, we have % G Hq{^^^). By 
appealing to the definition of gl we find that (c/. p.9p ) 

= / Vgl : \/(dx = / \/{gl + rj,) : VCd.T for aU C e H^i^")- (3-27) 

Hence, rj^ is the unique solution of p.27p in i7g(ri2). Since (Vg^) is bounded, taking ( = j/^ in 
(|3.27p shows that the sequence (77^) is bounded in i7g(ri^). Therefore, by Proposition IA.5I we get 

= 0{e) and, up to a subsequence, Vr;^ Vj,?7o in L^{QxY) for some 770 S L^(f2; i/p (^o))- 
By the periodic unfolding of the integral in the right- hand-side of (|3.27p (see Proposition IA.4P and 
a rearrangement of the terms in the resulting equality we get 

/ |V7?,|2 dx= ff U,{Vr,,) : UeiVri,) dxdy^ - (( U^iVgl) : U,{Vr^,) dxdy. 

Since the sequence (V5I) (resp. (Vr^e)) strongly (resp. weakly) two-scale converges, passing to the 
limit as e in the right-hand side of (|3.3p yields 

lim / |V7?ep dx- = - / / {Vg^{x) +Vyij{x,y)) ■.Vy7]oix,y)dydx. 
^^"Jn JJnxY 

Since 770 S L'^{^}; H^{Yo)), Vg{x) is constant with respect to y, and because tp satisfies p.9p (resp. 
p.Sp ). the right-hand side of p.3p vanishes, which means that Vr^e strongly converges to zero in 
L^(fi), as claimed. □ 
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Proof of Lemma \3.14\ Recall that we denote by C a positive generic constant that can be chosen 

independent of e. Since dYo is a Lipschitz boundary, there exists a family (E^) of continuous 
extension operators from W^'P{nl) to W^'P{ft) such that 

yv&W^'P{^l)-- l|V(ii;e77)lki.p(o) < C||V?7lUi-(ni) (3.28) 

(see e.g. [71 Appendix B]). Set r/* := E^{r]g\Qi) and := rjg — ij*. By construction we have £ 

Hq{Q^), and by assumption p.2ip and property p.28p we have Vrj* — > in LP{il,). Furthermore, 
by (|372^ we have 



/ Vwe : VCdx - / Vi]; : VCdx for all C & W^J'°°(f^°). 



(3.29) 



Thus, is the unique solution in Wo'^(rJ°) to a non-homogeneous Dirichlet problem with the right- 
hand side —At;* understood in the sense of distributions. By elliptic regularity we expect the a 
priori estimate 

||we||wi.p(oo) < C'||V7?*||lp(oo) (3.30) 

to hold. Since ||V77*| 1^,2(00) ^ this would complete the proof. For the precise argument we rescale 
Ue to domains independent of e and consider the components of the vector- valued map separately. 
Therefore, let i G {1, . . . , d} be arbitrary and for all ^ e Z'' the maps 



e-\U,vm,y) ife(e + >^o)cr!0 

else, 

Uu,y{ifJ){£,,y) if£(^ + yo)cl}^ 

1 else. 



where and (r?*)' denote the ith component of the corresponding vector. It is straightforward to 
check that 

Furthermore, by p.29p for all € the function is the unique weak solution in WQ^^iYt^) to 
the non-homogeneous Dirichlet problem 

-Alf = V-F« inFo, 
= on 9^0- 

By elliptic regularity we get the a priori estimate 

\\VHw-..(Y,) < qi^ilUnro) for aU C G 
as is shown in [25l Theorem 0.5] . By virtue of p.3ip this implies p.30p , and the proof is complete. □ 

3.4 Proof of Theorem \3A\ 

Recall that {g^,gl) ip and set (pi := gl and (^5° eg"; note that the deformation associated 
to If takes the form 

u^\d + ip = \d + ipl + fl=\d + egl + e^g]. 



Let us denote by the extension of wLi to as given in Proposition [521 We have 
||dist(V2„S0(n))||^.(f,) < C||dist(Vw,SO(ci))||^.(,,,) . 
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By assumption we have u = id on F C dfl C dfll, hence = id on F, and so the geometric rigidity 
estimate from Proposition 15.41 yields 



||Vu, < C||dist(Vu,SO(d))||^.(f,i) . 

By construction it is e'^gl = = — id, thus 

l|V^.!||i.(o) < C||dist(V^.,SO(d))||i.(^,). (3.32) 

We proceed now with ip'^, and start with the pointwise incquahty 

dist(V(u-</j^),SO(d)) < dist(Vu,SO(d)) + |V(^^|, (3.33) 

which is obtained as fohows. Fix x Cz Q and choose rotations Ri and i?2 in such a way that 
dist(V(u- (pi),SO(d)) = \V{u~ipl) ~ Ri\ and dist(Vu, SO(d)) = |Vit- i?2|. Then by the triangle 
inequality one has 

|V(u-(^i)-i?i| < |V(u-^i)-i?2| < \Vu-R2\ + \Vipll 

which is exactly p.33p . From Proposition 15.41 (recall that u ~ ipl = id on the boundary of 17°), 
applied via a scaling argument as in Proposition 15. 2[ and by the estimates (j3.33[) we get 

||V^°|U2(oo) = \\W{u - pI) - I\\L2^no) < C\\ dist(V(^. - pI), SO(ci))|U.(no) 
< C\\ dist(V«, SO(d))|U2(oo) + C||V<^i|U2(oo). 

Then via (|3.32p we finally conclude that 

l|V^°|li.(oo) < C||dist(V^.,SO(d))||i.(^). (3.34) 

Since (^^ = on F, we get 

MWrnin) < C||dist(V7.,S0(d))||i.(j,,) (3.35) 

by appealing to (|3.32p and to the standard Poincare inequality. By construction the displacement 
vanishes on dfl'^. Thus, (|3.34p and the high-contrast Poincare inequality (see Lemma [3. 61 ) yield 

MWmn) + lkV^?lli=(o) < ce^\\dist{\7u,S0{d))\\l.^^y (3.36) 

By assumption the stored energy functions Wi, i — 0, 1, are of class We- In particular, they are 
non-degenerate with constant 9 > 0; thus, we have 

II dist(V^., SO(d))||i.(^) + £-27|| dist(V^., SO(d))||i.(^,) 
< ||dist(Vu,SO(d))||i.(f,o) -HCe-2''||dist(V«,S0(d))||i.(f,i) < CA,{^). 

Now, the statement of the theorem follows from p.35p and p.36p . □ 

4 Asymptotics and F-convergence 

In this section we study the asymptotic properties of the total energy p.7p along sequences {(pe) C 
Hp ^{il) that satisfy the low-energy condition p.l4p . and for loads (L) C i^(il)'^ \ {0} that are 
relatively small in the sense of the condition (|1.2p . For reader's convenience we recall that this 
means that 

limsupAe<cxD and \imsnp ^ ( A^iip^) — [ l^-ip^dx] < oo, (4.37) 
e^O 6^0 H \ Jn J 

where := e^\\le\\L^{Q.) + £||^e||L2(no). 
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Let (g^^gl) be the e-splitting associated to ^ips via Definition 13. II Using the a priori estimate 
of Theorem 13.41 and Corollary 13.51 we find that 

^Pe = Ki^g'^ + s'^gl) and g^,ey g1,y g1 are equi-bounded in L^i^l), 

r o (4-38) 

Ae(v?e) = 0{\l), A,i^,) - dx = 0{\l). 

Jn 

We obtain several pieces of information from (|4.38p . First, clearly 

\We\\mm = 0{K), \W,\\mn) = 0(A,£"""^^'i>); 

and we immediately see that any low-energy sequence weakly converges to zero in H^(Vl) (since 
7 > by assumption). Furthermore, (|4.38p indicates that different kinds of behaviour are to be 
expected for ~ 1 and A^ — > 0, respectively. 

In the regime A,. — the infimum of the total energy converges to zero. Moreover, the sequence 
of displacement gradients iV'Pe) converges strongly to zero and thus the associated nonlinear strains 
■\/(/ + V((9e)^(/ + Viy9e) get arbitrarily close to the identity. Broadly speaking, this allows to replace 
the energy contributions of the soft and stiff component by their linearisations, i.e. 

I Wo{I + Vip,)dx^ Xl f Qo(eVg°)dx, 

£-2^/ Wi{I + V^,)dx^Xl [ Qi{Vgl)dx, (4.39) 

where Qi is the quadratic term in the expansion of Wi at the identity, i = 0, 1. Motivated by this, 
in Section 14.21 we study an appropriately scaled version of the energy, which leads to non-trivial 
limiting behaviour described by a quadratic energy functional. 

In the regime Ag ~ 1 the infimum of the total energy is of order one and large strains might 
emerge in the soft component, thus a full linearisation cannot be expected in this case. However, 
we argue in Section |4?3] that in the stiff component the nonlinear strain is still close to identity; and 
thus allows for a partial linearisation, in the sense that (|4.39p remains valid. 

4.1 Convergence of the load term 

In this paragraph we establish some convergence properties of the load term. Let us first refine our 
assumptions on the applied load. In order to exclude concentrations we suppose that there exist 
limiting loads l'^ £ L'^{D,) and /" G L'^{flxY) such that 



ll := XTh'^L f weakly in L^(n), 

2 (4-40) 
:= Ag ^eXe'e ^ ^^{XiV) strongly two-scale 



Remark 4.1. // the first of the assumptions (|4.37p is satisfied, then the components X'^^e'^l^ and 
Xj^eXeh o.Te equi-bounded in i^(f2). Thus, only the assumption that the latter component of the 
load is strongly (rather than just weakly) two-scale convergent introduces an additional constraint. 
Conversely, if a sequence (l^) C L'^{il) satisfies (j4.40p with respect to a given bounded sequence (X^), 

then e'^||Ze||L2(o) -|- £||/e||L2(oO) Ae. 

It is convenient to introduce the following notation for the load term: for {g^,gl) e Ge{^) and 
(ff°,5^) G Gr,o we set 

C%gl gl) := / ?° • .g° Ax + ( l] ■ g] dx. 



C'ig',g'):= I (I f-g'dy + l^-g^] 
Jn \Jyo / 



(4.41) 
dx. 



Clearly, we have 

^Jje- dx = C^g^gl) for {g^gl) ^ Ar Vs 
Assumption (|4.40p is tailor-made for passing to the limit in 
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Lemma 4.2. Let C L'^{n) \ {0} satisfy tf4A0\ . Let ig°,gl) G Gei^) and suppose that 

eVgl A v,gO, Vgl - Vg^ with {g\g^) e Gr.o- Then 

lim£-(77°,r;i) = C'{g\g'). (4.42) 

Proof. Since g° = on 9ri° (resp. = on F), the convergence of eVg^ (resp. V^^) implies 
that g^ g^ weakly two-scale and gl g^ strongly in L^. Thus, (|4.42p follows by (|4.40p and 
Proposition IA.41 

□ 

4.2 Convergence in the small-strain regime (Ag — )■ 0) 

As is mentioned above, for — > the asymptotic behaviour is trivial: the infima of p.7p and the 
associated low-energy sequence [ip^) converge to zero as e — s- 0. In particular, this applies to the 
trivial displacement — 0. For a more accurate description of the asymptotic behaviour we study 
the scaled elastic energy \~^K^. We prove that it F-converges to the limit energy 

^smaii^^o^^i) j^^ (^J^Qo{yyg')dy + QT'i^g')^ dx, (g°,.gi) e L\n;H^{Yo)) x Hl„{n) 

with respect to the weak two-scale splitting convergence introduced in Definition 13.81 Recall that 
Qi is the quadratic form in the expansion of Wi at the identity (see Definition 12. ip . We denote by 
Qhom . jgdxd [0, od) thc rcsult of "perforated" homogenisation for Qi which is determined by the 
formula 

g^°'"(^):= niin / Q,{F + WHy)) dy. (4.43) 

Theorem 4.3. Suppose that > and A^ — > as e — > 0. Then the sequence A^^A^ converges to 
^smaii following sense: 

1. (Lower hound). Let {ipe) C q{^) and {g'^,g^) G Gr,a- The following implication holds: 

A7Ve^-^'(5°,ff') liminfA-2A,(^,)>A™-"(5",.gi). 

e— >0 

2. (Existence of recovery sequences). Let {g^,g^) G Gr,o- There exists a sequence (ip^) C Hp Q{il) 
and ijj S Gr.o such that 

AJ Ve 'A' (/, g\i^) and Imi Xj'Uv^e) = A™^"(.gO, g^). 

Before proceeding to the proof, we use the convergence result above to gain a refined asymptotic 
analysis of the minimisation problem for p.7p . Namely, we show that the asymptotic low-energy 
behaviour of p.7p is effectively described by the functional 

Let {9^,gl) € C/r.o(^^xy) be the unique minimiser of the above limit energy functional and let -0* 
satisfy 

/ g'r°(V3i(a;))d.T= // Q,hgl{x)+VyM^,y))dydx. 
Jn J JnxYi ^ ' 

Proposition 4.4. Suppose that A > 0, A^ 0, and that (l^) G L'^{il) \ {0} converges to and l^ 
in the sense of (|4.40p . 

1. (Convergence of infima). The convergence 

hm inf ^(a,{p)- f k-pdx) =A^'^^"i9lgl)-C'{gl9l). (4.44) 

holds. 
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2. Let {ips) C Hy o(i^) be a sequence satisfying 

A,{^,)~ [ l,-ip,dx< inf (A,{^)-[l,-^dx]+o{\l). (4.45) 
Jn veH^o(n) V Jn J 

Then: 

(a) (Weak two-scale convergence of minimising sequences). The convergence 
holds. 

(b) (Strong convergence of minimising sequences). Let {(p^) be a recovery sequence, i.e. 
Xj^ipe '-^^ i9*,9l,'4'*), and consider {g°,gl) ^ K^fe, ide^gl) K^Ve, then 

\\9e - 9e\\Lp{no) + ||£V5° - e\/g°\\LP{no) + \\gl ~ ~g\\\w^.Hn.) -> 
for alll <p <2. 
Wc prove the proposition in Subsection 14.2.21 

4.2.1 Proof of Theorem llTSl 

The following proof is a slight variation of an approach presented in [31j and combines convex 
homogenisation methods with a linearisation argument inspired by some of the methods in [25] . 
The core of this approach is the following statement (c/. [3TJ Theorem 5.2.1]). 

Lemma 4.5. Let W be a stored- energy function of class Wg and denote by Q the associated quadratic 
form in the expansion ofW at the identity. Let (i^^) C L'^{n), F € L^{flxY) and let Cg be a sequence 
of positive numbers converging to zero. 

1. LfF.^F, then 

liminf cjM VF(/ + CeFe)dx > / / Q{F{x,y)) dx dy. (4.46) 
-^^0 Jn JJnxY 

2 

2. If Fi; — > F and limsupg_^Q C£||i^e||/,oc(-f2) = 0, then 

limc^M VK(/ + CeFe)d.T = // Q{F{x,y))dxdy. (4.47) 
Jn JJnxY 

Proof. The proof relics on the weak lower semicontinuity (resp. strong continuity) of convex integral 
functionals and the property that W possesses quadratic behaviour near the identity. Therefore, let 
us recall the assumption (|2.6p : 

W{I + G)=Q{G)+r{G), where r(G) = o(|G|^) as |G| ^ 0. (4.48) 

Without loss of generality we assume that for all matrices G we have 

\r{G)\<p{\G\)\G\\ 

where p : [0, +(X)) — > [0, oo] is an increasing function with limg^o p(<7) = p(0) = 0. To shorten the 
notation we define for G E L^{^1) the functionals 

/e(G) ■.= c~^ [ WiI + c,G{x))dx. 
Jn 
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Taking advantage of the property of the unfolding operator lA^ referred to as "integral identity" in 
Proposition |A21 we find that 



L{G) = c-^ ff W(l + cMe (G) (x, 2/)) dx dy . 

JJr''xy ^ ^ 



(4.49) 



Recall that the support of the periodic unfolding U^[G) associated to a map G G is slightly 

larger than fixF. Wc denote by D a bounded set such that supping (G) C fixF for small 
enough e and all G G Furthermore, wc tacitly extend maps defined on OxF to M'^xF by 

zero. 

Proof of (|4.46p . In order to apply the quadratic expansion (|4.48p we "truncate the peaks" of and 
consider 



Fe{.x,y) := 



\U,{F,){x,y) if U,{F,){x,y) 
[ otherwise. 
We claim that this modification docs not affect the limit, i.e. 

F weakly in {M.'^xY) 



< c] 



-1/2 



(4.50) 



which can be seen as follows. By construction the sequence (F^) is bounded in L'^{Q,xY) and is 
therefore weakly relatively two-scale compact (see Proposition IA.4P . Thus, it suffices to show that 
Fe — Ue{F^) in L^{tl'KY). Indeed, by the Holder inequality we get 



U,{F,)~F, 



< 



U,{F,)~F, 



,1/2 



(4.51) 



where 



{F.^WSF,)] = {{x,y) : \WSF,)\>c~''^] 



Now, the Chcbyshev inequality yields < Ce\\F^\\'j^2i^Qy Thus, the right-hand side of (I4.51|) goes to 
zero as £ — !• and (|4.50|) follows. 
Next, we argue that 



liminf /^(Fg) > liminf 



£^0 



6^0 



Q{Feix,y))dxdy. 



(4.52) 



'xY 



By assumption, W is minimised at the identity. Since Fg{x, y) equals either Ue{Fi;){x, y) or zero, we 
get W[I + CeUe{Fe){x, y)) > W[l + CeF^ix, y)) . With (|4.48l) this yields the pointwisc lower bound 

c-^w(l + cMeiFe){x, y)) > Q{F,{x, y)) - p{c,\F,ix, y)\) \F,{x, y)\\ 

We integrate both sides of the last inequality over R'^xy. In view of (|4.49p . the left-hand side of 
the result equals Ie{Fe). Since p(ce|Fe(a;, y)|) < p{^/c^) — o(l) and because the sequence (Fg) is 
bounded in L^, we get 



lim 



pU\F,{x,y)\]\F,{x,y)\^ Axdy = 0, 



and (|T5^ follows. 

Now, (|4.46p follows from (|4.52p . (|4.50p . and the lower semicontinuity of convex integral functionals 
with respect to weak convergence (see e.g. [H]). 

Proof of (|4.47p . By appealing to the expansion (|4.48p we obtain the pointwise estimate 
c-^w(l + cMe{Fe){x,y)) <Q{u,{F,){x,y)) + p{c,\\F,\\L,.(n))\Ue{Fe){x,y) 
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Since U^{F^) is bounded in (as a two-scale convergent sequence) and because of the assumption 
limsupg Ce||i^e||i=o(Q-) = 0, the pointwise inequahty impUes that 

hnisup/e(Fe) < Hmsup / / Q{Ue{F^){x,y))dxAy. (4.53) 

By assumption U^{F^) strongly converges in i^(R''xF) to F (extended by zero to M'^xF), hence 
the right-hand side of (|4.53p equals JJ^^y Q^^) da-'dy and we get 

limsup4(^e)< // Q{F{x,y))dxdy. (4.54) 

Now (|4.47l) follows from (|4.54l) and the lim inf-incquality (I4.46|) . □ 



Proof of Theorem 14.31 (Lower bound). Wc only need to consider the case when '^A^{ipg) 
converges to a finite number as e — > 0. By appealing to Proposition 13.101 we pass to a subsequence 
(not relabelled) such that A^Ve (5*',3\ V') with {g",g^,i;) e GrfiinxY). Set 

F^> := X-'x^e^'Pe, Fl X-^e'^xl'^Ve- (4.55) 
By definition of the e-splitting, we have 

Ff' = eVg'i + e-'x^gl , Fl = gl , 

where {g^,gl) «^ X~^ip^. Since the differences xli^) ~ X^ix/e), i = 0, 1, converge to zero in L^(57) 
as e — >■ 0, we get via Proposition lA. 41 



,g°(x,y), Fl ^ x\y){^9\^) + ^y4'{x,y)). 
It is easy to see that 

A72A,((p^) = A7^ / Wo{I + \eF°{x))dx + \-'^e~'^'' [ Wi{I + X^e^ Fl{x)) dx. (4.56) 



Indeed, by assumption Wi vanishes at the identity, therefore Wi{I + G)x = + xG), i = 0, 1, for 
all matrices G and x € {0, !}• Now (|4.56p follows from the definition of {F^,Fl) as a decomposition 

of (fie- 

Finally, we apply the lim inf inequality from Lemma 14.51 to each of the two integrals in the 
right-hand side of ()4.56p and get 

liminf A-2a^(^^) > ff Qo{Vyg°{x,y))dxdy+ [[ Q^hg^x) +Vyi:{x,y)) dxdy. 

JJnxY J JnxYi ^ ' 

Since the second integral is bounded below by Q5'°'"(V(7^(x)) dx, the lower-bound statement of 
Thcorcml4. 31 follows. 



Proof of Theorem l4.3l, (Existence of recovery sequences). Let V' be a minimiser in iP'i^; H^^IY)) 
of the functional 



nxYi 



(v5i(x)+Vy7/.(x,y))dxd?/. (4.57) 



Since the domain of integration in (j4.57p is flxYi, wc can define ip in ^IxYq at will. In particular, we 
assume that -0 G Qr.oi^xY). Let (ip^) C i?pQ(ri) denote the sequence associated to {g''\g^,ip) via 
Proposition 13.131 with ~ A^ and set ipi, = Xe^e- As in the proof of the lower bound, wc consider 
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the decomposition {F^,F^) of the displacement gradient Wife according to the formulae (|4.55p . so 
the identity (|436| holds. By Proposition EH] (c/. ([3:23| . (fSTM]) ) we have 

lim(A,||i^O|Uoo(n) +A,£^||^^,i|U^(o)) =0. 

Now the second statement of Lemma 14. 5[ applied first with ^ X^, F^ — and then with 
c, = A^e'^, Fe = F^, yields 

lim Aj2A,(^e) - // Qo{Vyg°ix,y))dxdy + [[ Qihg\x) +Vy^{x,y)) dxdy. 

JJnxY JJnxYi ^ ' 

This concludes the proof, in view of the fact that the second integral in the last formula equals 
LQ^"(V5)dx. □ 

4.2.2 Proof of Proposition [4^41 

The proof is divided into seven steps. In Step 1 we deduce convergence of minima and weak two-scale 
convergence of minimisers by appealing to standard arguments from F-convergence. The goal of the 
remaining steps is to improve the convergence to strong two-scale convergence. In a nutshell, the 
idea is to "gain strong convergence by strict convexity" . In order to apply this paradigm we have 
to replace the non-convex energy by its linearised version, the main part of the related argument 
is contained in Steps 3-5. In Step 3 we show that for minimising sequences the linearised energy 
converges as well, at least after "truncating peaks". To handle the truncation error, in Step 4 we 
modify the minimising sequence on small sets to gain equi-integrability. Finally, in Steps 6 and 7 we 
compare the minimising sequence with the recovery sequence by estimating the differences between 
the original, the modified and the recovery sequence. 
For convenience we introduce the abbreviations 

eo := A--"(.9°,gi), h ■= C'^iglgl), mo eg - ^o- 
Furthermore, let {g'^,gl) denote the e-splitting associated to A~^(^e. 

Step 1 (Convergence of the energy and weak convergence of minimisers). As well as for classical 
F-convergence, the "lower bound" and "recovery sequence" properties in Theorem 14.31 together with 
the compactness result in Proposition 13.101 implv the convergence of the infima in (|4.44p and of any 
minimising sequence to the minimum and the minimiser, respectively, of the limiting functional. 
Indeed, since Ag{ip) — J^l;. ■ ipdx = ioi (p = 0, one has inf^^g^i ^(j^-j (Ae(iy9) — J^^ 1^ ■ ipdx^ < 0. 
Then the condition (|4.45|) implies via Corollarv l3.5l that limsup£_j.o X~'^Ag{(pe) < oo. By Proposition 
13. 101 there exists a triple p^, V') G Gr o{^xY) such that, up to a subsequence, one has X~^(ps 

By the assumption (|4.45p . Theorem 14.31 (Lower bound) and Lemma 14.21 we have 
liminfj_^o infi/^eifi ^(O) {Aei'p) ~ ■ fdx) > mo. In turn, the "recovery sequence" property 
of Theorem 14.31 and Lemma l42l vicld the opposite inequality, and (|4.44p follows. By the uniqueness 
of the minimiser we get {g'^,g^,ip) = (<?* , 5*, "0*)! ^-^d hence 

, _ 1 sp2 ,0 1 / \ 

K (g*,9*,v*)- 

Step 2 ((Lower semi-) continuity of quantities associated to the linearised energy). In this step we 
gather some continuity properties for the linearised energy, the associated bilinear functional and 
the load term. For matrix-valued functions G'°,F°,G^,F^ e L'^{Vl) define 

Q'{G\G^)~ f Oo(G'°)dx+ / QiiG^)dx, 

B'{G'^,G';F'^,F') := f Bo(G°;f^°)dx-t- / Bi{G\F^)dx 
Jn° Jni 
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where Bi, i = 0,1, denotes the bihnear form associated to the quadratic form Qi, i.e. 

1 



yA,BeR'''"' : B,{A-B) := -{Q,{A + B) - Q,{A) - Q,{B)) . 
The following (lower semi-) continuity properties hold: If 

G^-GO, Gl^G^ and F° A F° , A F\ 

in L^, then 

hminf Q-(G°,Gi) > (^^^ Qo(G°(x, y)) dy + ^ Q,{G\x,y)) dy'j dx. 



(4.58) 



>0 



n \JYo 



\imQ%F^,F^)^ I (I Qo{F°{x,y))dy + J^Q,{F\x,y))dy'jdx, (4.59) 



\miB%GlGl;FlF})= [ ( [ B^{G\x,y),F\x,y))dy (4.60) 

JS2 \JYo 

Bi{G\x,y),F\x,y))dy\ dx. 



The statements (|4.58p and (j4.59p are well known and follow by "lower semicontinuity with respect 
to weak two-scale convergence" of convex integral functionals and "continuity with respect to strong 
two-scale convergence" of continuous integral functionals (see e.g. [37|). The statement (|4.60p is due 
to the fact that it contains only products of strongly and weakly two-scale convergent quantities. 

Step 3 (Convergence of the truncated linearised energy). We claim that 

lim ( / Qo{eVg°)dx+ [ Qi{Vgl)dx - C'{g°,gl)] = mo (4.61) 

for a sequence {X^) of measurable subsets of il with 0. 

For purposes of our argument let consist of those x E fl where min{|eV(7°(a;)|, |V(7£(a;)|} > 
AJ^^^- Recall that the sequences (eVg^) and (V.g^) are bounded in L^(0). Hence, in the same way 
as in the proof of Lemma 14.51 one has 

|X,|->0, %Vg°^V^.g2, X^Vg^ ^ VgJ + V,V*, (4.62) 

where Xe denotes the indicator function of f2 \ X^. Due to the quadratic expansion for Wi, i = 0, 1 
(see ()2.6p ) and the fact that e'^Vg^ strongly in L^, we find, in the same way as in the proof of 
Lemma l475l that 

1 

a! 



Mfe)- f le-Vedx) >{[ QoieV g'o) dx + [ QiiVg^dx] -C'{glgl)+r, 

(4.63) 



where — > 0. Since the left-hand side of (|4.63p converges to mo, we get 



mo > lim sup 

£^0 



Qo{e^gl)dx+ f QiiVgDdx] -C'iglgl) \ 
Jnl\x, J J 

Qo(eVgg)da;+ / Q,{Vgl) dx] ~ C%gl gl) \ 
Jni\x, J J 

[ ( [ Qo(V,.g") dy + [ Qi{Vgl + V,^,) dzA dx - C\g\g^), 

JQ. \JYo JYi / 



> lim inf 



> 
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where we used (|4.58p . (|4.62p and (|4.42p in order to pass to the hmit. Since the last expression in 
the previous estimate equals toq, the assertion follows. 

Step 4 (Equi-integrable decomposition) . In order to motivate the following steps, let us remark that if 
the linearised energy Q^(eV(7°, Vg^) — [g^^g].) converged to the minima mg, then we could obtain 
strong two-scale convergence of eS/ g^ and Vg^ by appealing to the paradigm of "strong convergence 
by strict convexity". However, (j4.42p ignores the energy contribution emerging in the subset X^. 
Clearly, if we knew that 

lkV<?°||i2(nonx,) + W^glWh^ninx,) ^ 0, (4.64) 

then (|4.42p would also hold for = 0. Unfortunately, (|4.64p is not valid for general sequences. To 
overcome this difficulty we modify (g^^gl) on small sets in such a way that the resulting sequence 
is equi-integrable and (|4.64p is satisfied. To make this precise, we claim that for a subsequence 
(not relabeled) there exist sequences (g°) C H^{n), {§1) C ^{Q) such that gl satisfy (|3.8p . the 
convergence statements 

eVg° - eV5° ^ 0, \7gl - Vgl ^ 0, (4.65) 

||eVg° -£V5°||LP(n) ^ 0, \\Vgl - VglWL.in) ^ (4.66) 

hold, and 

lim (Q^'(£V5°, Vgl) - C'{gl gl)) = mo. (4.67) 

The argument relics on the well-known Decomposition Lemma by Fonscca, Miillcr and Pedregal 
(see [22]). By appealing to Lemma [4.61 and passing to a subsequence (not relabelled), we find that 
there exist sequences (g°) C 77o(ri), {gl) G q{^) and "exceptional" sets 

X,:^{x^VL: g^,{x) ^ f,{x) or gl{x) ^ gl{x) } 

such that the sequences (|eVg°p) and (|Vg^p) are equi-integrable and \X^\ — > 0. Clearly, by 
construction, (|4.65p and (|4.66p hold whenever 1 < p < 2. It remains to prove the key identity (|4.67p . 
By virtue of (j4.65p , the lower semicontinuity property (|4.58p and Lemma 14.21 it suffices to argue 
that 

lim sup {Q'ieVglVgl) ~ C'{t,gl)) < ^0- 

Let Xe and Xe denote the characteristic function of (the set from step 3) and X^, respectively. 
To simplify the notation we use the shorthand := (eVg^, Vg^) and := (eVg^, Vg^). Taking 
into account (|4.6ip . it suffices to show that 

limsup (Q'iGe) - Q'iH - Xe)Ge)) < 0. 

By the quadratic property of and by the (pointwise) orthogonality of Xe and 1 — Xei fii^d 
Q'{G,) - Q^((l - Xe)G,) = Q'ixeG,) + Q^((l - Xe)G,) 

- Q'(Xe(l - Xe)Ge) - Q"((l " Xe)(l ^ Xe)G,). (4.68) 

By definition of the set X^ we have (1 — Xe)G£ = (1 — Xe)Ge, thus we can simplify the sum of the 
second and fourth term in the right-hand side of the equation (|4.68p as follows: 

Q^((l - Xe)Ge) - 2^(1 - Xe)(l " Xe)G ,) = Q^(Xe(l - X.)G,). 

Using the fact that the third term in the right-hand side of ()4.68p is negative, we arrive at the 
inequality 

Q'{Ge) - Q"((l - Xe)G,) < Q'iXeG,) + Q^(Xe(l - Xe)G,). 

By construction the sequences (xe(l — Xe)) and (xe) converge to zero boundedly in measure as 
e — > 0, thus, by equi-integrability the right-hand side of (|4.68p vanishes in the limit as e — and 
the argument is complete. 
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Step 5 (Comparison with the recovery sequence). Recall that denotes a recovery sequence associ- 
ated to {g'^,gl,'4'*) and ^e^Ve- Let ig^,gl) denote the modification of {g^,gl) according 
to (|4.65P " (|4.67p . In order to monitor the difference between the modified sequence {g^,g\) and the 
recovery sequence [cf^^g^), we define 



t---t~~gl vl^^gl^gl- 



Wc claim that, up to a subsequence, 



/ |symV(e^°)pdx+ / jsymV^^pda; 
In order to show (|4.69p . we set 



^ as £ ^ 0. (4.69) 



eo := lim Q^(£V5^^,V5,^). 

e— >0 

Note that the above limit exists since {gg,gl) is a recovery sequence. Furthermore, by (|4.67p and 
Lemma [4?2l we deduce that 

lim Q^eVgl Vgl) = lim QHeV~gl V~gl) = eo- 

e— >0 6— >0 

Clearly 

Q^e\/7t,\/fjl) = Q^e\/g°,\7gl) - 2B%e\7f,\/gl;e\/g°,\7gl) + {eV f,,^ ~gl) , 
where we can pass to the limit in the bilinear term by ()4.60p : 



ix) + VV'*(a;, y)) dy I dx = eo 



B'ieVf^Vgl-e^f.^V-gl)'^ I ([ Q^{V ygl{x,y)) dy + ( Qi(Vgi 

Jn \Jyo Jyi 

Therefore, Q^(eV^°, Vfjl) and since the quadratic forms Qq and Qi are positive definite on the 
space of symmetric matrices, (|4.69|) follows. 

Step 6 (Convergence of the stiff part). In this step we prove that 

\\gl-gl\\w^.P(n)^Q as e ^ 0. (4.70) 

We set nil := gl — gl, and assume without loss of generality that 3/2<p<2. By construction we 
then have J^o "^vl ■ V^dx = for all ( € i?o(f2g) which allows one to apply Lemma r3.141 once it is 
shown that 

||?7£ llwi-p(ni) 0- 
Note that rjl ^ gl — gl + fjl and consequently 



\\vl\\w^-p{ni) < \\gl -5ellwi.p(0) + \\vl\\w^-H(il)- 

The first term in the right-hand side of the last inequality vanishes up to a subsequence due to 
(|4.66p and the property that gl = gl on F. Regarding the second term, it suffices to show that 
Il^ell//i(f2i) — i> as e — >■ 0. To this end, let fjl denote the extension of ^^Ini to il according to 
Lemma [3T5] Using the boundary condition fjl = fjl = on F, the Korn inequality and p.l2p we get 

WvlWH^n) < c'||symV7)^^||L2(a) < c"|| sym V^^^||l2(oi). (4.71) 

Now, the left-hand side of (|4.7ip controls ||f^^||iji(oi), while its right-hand side vanishes by (|4.69p 
up to a subsequence. Since the above argument applies to each subsequence, (j4.70p holds for the 
entire sequence. 

Step 7 (Convergence of the soft part). In this step we prove that 

ll5°-5elU^(no) + ||eV(g°-g°)|Up(no)->0 as e -> 0. (4.72) 
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The argument is similar to that of Step 6. We set 77° := — 3°, where is defined according to 
(|4.65p and ()4.66p . Since 77° = on dfl^ we have, by applying the Poincare on Yq and rescaling: 

||'7ellLf(no) < c||esymV?7°||LP(no). 

The identity 7]^ = .g^ - .g° + 77° yields 

WVeWLPino) < c' \\esymWT]°\\ LP (no) < c" (HesyrnVg" - esymWg°\\Lp{no) + ||esym V77°||l2(qo)) . 

The right-hand side of the last inequality vanishes up to a subsequence due to (|4.66p and (|4.69p . 
Since this argument applies to each subsequence, ()4.72p holds for the entire sequence. 

□ 

Lemma 4.6 (Decomposition Lemma, see [35], ^7\). Let v G Hr '^"'^ (''^e) ^ sequence with 
V in H^{fl). Then there exists a sequence (ve) C H^(n) such that the following properties hold 
for a subsequence of (not relabelled): 

1. V, -^v in H^{n); 

2. = in a neighbourhood of dfl; 

3. (|Vwep) is equi-integrable; 

4. |{ 2; G : Ve{x) ^ Ve{x) }\ 0. 

4.3 The finite-strain regime (A^ ~ 1) 

In this section we consider the case ^ 1, when, unlike in the small-load regime, the strain emerging 
in the soft component is not necessarily small but only bounded as e — ^ 0. As a consequence, passing 
to the limit as e — >■ does not result in a linearisation of the functional on the soft component so 
that we have to deal with the full non-convex problem. It is well-known from the work by S. Miiller 
|30| and A. Braidcs [7], where the homogcnisation of non-convex integral functionals is discussed, 
that minimisers of periodic non-convex integral functionals tend to oscillate on scales larger than 
the period of the material. As a consequence of this, the homogenised energy density is given by 
a mutli-cell homogcnisation formula that involves minimisation over an infinitely large ensemble of 
periodicity cells, in contrast to the one-cell homogcnisation formula of the convex case. For the 
same reason a lower scmicontinuity result for non-convex integral functionals with respect to weak 
two-scale splitting convergence docs not hold in general. 

To overcome this difficulty we prove F-convergence of with respect to the stronger two- 
scale splitting convergence as defined in Definition 13.81 which forces the gradients V^g" to oscillate 
essentially on the scale e. However, wc find, as an associated drawback, that in general the sequence 
(Ae) is no longer equi-cocrcivc, and thus wc do not gain convergence of the related sequence of 
infima and almost minimisers by the standard machinery of F-convergence. While in general in 
the non-convex situation the strong two-scale convergence of minimisers cannot be expected, in our 
context we are able to prove convergence of the sequence of infima to the minimum of the F-limit by 
following the arguments developed in [M]. We also refer the reader to [M] for a detailed discussion 
of the related issues. 

4.3.1 Modified setting and main result 

Next we describe the problem setup in the finite-load regime. Since in this case no linearisation takes 
place on the soft component, we drop the assumption that Wq admits a quadratic expansion at the 
identity. However, we now assume that Wq satisfies a quadratic growth condition at infinity and 
is quasiconvex. Note that the combination of both conditions implies that Wo is locally Lipschitz 
continuous. More precisely, we assume that for some positive constants 9,Q > : 

• Wo is a quasiconvex map from R'^^'' to [0,oo) such that Wo{I) = 0; 
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• 6'dist^(i^,S0(d)) < Wq{F) < 9(1 + \F\^) for any F € W^""^; 

• \Wo{F + G) - Wq{F)\ <e(l + |i^| + |G|)|G| ioi any F,G eM.'^'"^. 

Under the above conditions we prove that the sequence (A^) F-convcrgcs to the hmit energy 
A'-8^(50,5i) (^J^WoiI + Vyg')dy + Q'r\Vg')^ dx, {g\g^) G L\n-Hl{Y^)) x Hl,{a), 



with respect to the strong two-scale convergence as introduced in Definition 13.81 Recall that Q\°^ 
is the homogenisation of the quadratic form from the expansion of Wi as defined in (14.43^ . Further- 
more, we analyse the total energy including the force term. 
Without loss of generality we assume that Aj = 1 

Theorem 4.7. Suppose that Wq satisfies the assumptions listed above and Wi is of class Wg. Then 
Ag converge to the functional A''^''*^" in the following sense: 

1. (Lower hound). Let {ipe) C Hy q{^) and {g'^,g^) G Qrfl{i^'xY). Then 

Ve^^ig^a') ^ limiufA,(^,)>A'-s<=(gO,5i). 

2. (Existence of recovery sequences). Let {g^,g^) G Qy.q{^X-Y). There exists a sequence {ife) C 
Hy i^(^) such that 

5^,2 (gO^^i) lim A,(^,) = A'-s<=(/,gi). 

Suppose that [l^) C L'^{il)'^ converges in the sense of (|4.40p . 

3. (Convergence of the load term). Let (ife) C q{^1) with limsup^ A£((y9e) < oo. Then 

{g",g') ^ lim / ^,-kdx = £"{g",g'). 

4-. (Convergence of infima). The convergence 

lim inf (kJ^,)~[^,-l,dx^= min (A^^^^c^ o 1n o l^^ 

holds. 

Remark 4.8. By quasiconvexity of Wq and the argument of \14\ Section 9], the limit functional 

admits a minimiser, which is a property that is used in the argument for (Convergence of infima). It 
should be mentioned that the non- degeneracy condition is stable under quasiconvexification. Indeed, 
Zhang proved in fSSf that the quasiconvexification of the function "^(F) ~ dist^{F,SO{d)) is bounded 
below by c'^* itself where c' is a constant that only depends on the dimension d. For this reason the 
quasiconvexity assumption on Wq is not too restrictive. 

Remark 4.9. Our result can be generalised to non-quasiconvex Wq in a straightforward way by 
appealing to a lower semicontinuity result developed in \14^ . However, the restriction to quasiconvex 
Wq in Theorem \4-7\ makes the present paper self-contained. In a nutshell, in the authors show 
(for continuous W satisfying standard growth conditions) the lim inf inequality 



lim inf / W{sVve)dx> QW{\/yv)dydx 
J no JJnxYo 
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along sequences (v^) C H^{fl^) with Ve ^ v and sVve VyV where v S L"^ {U; Hq (Yq)) . Here QW 
stands for the quasiconvex envelope ofW. Further, they construct recovery sequences (converging in 
the same topology) such that 

lim / W^(eVue)da; // QW{\/yv)dydx. 

The key difference between the two approaches lies in the notion of convergence. In contrast to the 
convergence type used in Theorem \4. 7| which corresponds to strong two-scale convergence of e'Vvg, 
the results in fi/^f only require weak two-scale convergence of eVv^. In this way, the quantity eVv^ 
might oscillate on scales larger than e, which makes the analysis more subtle but still allows one to 
carry out the quasiconvexification. 



4.3.2 Auxiliary statements 

For the proof of the theorem and the subsequent discussion it is convenient to decompose and 
^large j^^^q ^^y^ energy contributions associated to the soft and stiff component. Namely, for {g^,gl) G 
H^in°) X H^^gin) and ig°,g^) S Gr.oii^xY) wc define 

A'M)--=f W,{I + eVg",ix))dx, Al{gl) :^ s-'-^ f Wo{I + e^Vgl{x)) dx, 

Aj](,9°):=// Wo{I + Wyg'{x,y))dydx, Aj(ffi) := / QT^V^'G^)) dx. 

JJnxYo Jn 

Clearly, we then have A^'^''s°{g'~' , g^) = Aq((;°)+AJ((7-'^). The following estimate is a direct consequence 
of the local Lipschitz continuity of Wq and the a priori estimate of Theorem [ 



Lemma 4.10. There exists a constant c > such that for any ip G Hp Q{il) and any e > the 

inequality 

A,{cp) - (AUgl) + AligD) I < c(l + A,i^))s'' 

holds, where {g^,g\) ip. 

Let us briefiy sketch the argument for the convergence proof for the infima in Theorem 14.71 To 
this end let be a displacement and (g'^,gl) the associated e-splitting. By construction we have 
(7° e i7Q(f2g) and by definition of the unfolding operator, the function Ug{g^) can be identified with 
a map in L^{fl; H^{Yo)). The proof of (Convergence of infima) relies on the following lemma. 

Lemma 4.11. Let ip E H^ ^{ft) and {g^,gl) «^ (p, then 

W_(gt) G L\^; Hi (Yo)), VyU,{g',) = £Z^e(Vg°), and A°(g°) = A° (uM)) ■ 

Together with the decomposition indicated in Lemma I4.10[ the above identities allow one to 
minimise the energy associated to the soft component at an intermediate stage {i.e. for e > 0). In 
this way we get for positive e > a sufficiently good lower bound for Ag{(pg), which is essentially of 
the form 

min A«(5")-// g' ■ l' dydx + Al{gl) - f gl ■ lldx. 

Importantly, only the energy contribution associated to the stiff component depends on e. Since 
the behaviour of Ag on the stiff component is nice, in the sense that the corresponding energy 
contribution A^ is lower semicontinuous with respect to weak two-scale convergence, we can pass to 
the lim inf and obtain the convergence of the infima. 
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4.3.3 Proof of Theorem l4?fl 

The proof is divided into five steps. We consider the behaviour of the elastic energy on the soft 
and stiff component separately. The argument for the stiff part, where the linearisation takes place, 
is similar to the one used in the proof of Theorem 14.31 For the soft part we make use of the 
strong two-scale continuity property of continuous integral functionals satisfying standard growth 
conditions (see for instance |37[). The crucial element in our proof is Step 5 below, where we 
prove the convergence of infima, by appealing to the equi-coercivity of with respect to the weak 
two-scale splitting convergence introduced in Definition [ 



Step 1 (Strong two-scale continuity of h.1). Let if^ [9^ id^) E^nd (g^,gl) (p^, then 

limAO(gO)= //' Wo{I + W.yg''{x,y))dydx. (4.73) 



In order to sec this, notice first that one has 



/ W^o(/ + £V5°(a;))dx= // Wo(l + U,{eVg°,)ix,y))dydx, 
Jn° JJnxYo ^ ' 

where Z/^e denotes the unfolding operator (see Definition lA.ip . Using the strong two-scale convergence 
of eV.g° and the continuity of Wp we have, for a subsequence, 

t¥o(/+Z^e(V^e)(a;,y)) ^ VKo(/ + Vy5°(x,2/)) a.e. in l^xFo- 
Furthermore, the quadratic growth of Wq at infinity yields the bound 

W^(l + eU,{\lgl)(x,y)^\< c\l + eU,{^ g%x,y)'\ 



for c > 0, where the expression in the right-hand side converges strongly in L^(i7xyo)- Our claim 
follows (first for a subsequence) by dominated convergence. Since the above reasoning applies to 
any subsequence of a subsequence, (j4.73p holds for the entire sequence. 

Step 2 (Proof of the lower hound). Clearly, we only need to consider the case when h^(ip^) con- 
verges to a finite number. Let (g^,gl) ip^. Since ip^ ^ (s'^iff^) by assumption, we find using 
Proposition 13.101 that 

eWg", 4 Vyg"{x,y), Vgl - Wg\x) + \/y^{x,y) 

for some -0 G L'^{^', H^{Y)) and up to a subsequence. 

By virtue of Lemma [4.101 and Step 1, it suffices to prove that 

liminf£-27 f Wi{I + e'''^gl{x))dx> [[ Q^""^ {S/ g\x)) dx . 
Jn JJn 



But this follows in the same way as in the proof of Theorem 14.31 

Step 3 (Proof of the existence of recovery sequences). Let tp £ QYfi{^y~Y) be a minimiser of the 
functional 

Qi(\Jg^{x)+\Iyil){x,y) \ dydx, 



IClxYi 

and {(f^) C Hy q{^) denote the sequence associated to (g'^jg^) and tp via Proposition 13.131 Then 
(fie ^ (g^Td^) Emd in particular we have 

eVg", 4 \/yg{x, y), Vgl 4 Vg{x) + Vyi,{x, y), 
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where {g'^,gl) cp^. By Lemma [4.101 it suffices to show that 

limAO(gO)=// W„{I + Wyg"{x,y))dydx, (4.74) 

limAi(gi)= JJ^Q\°"\Vg\x))dx. (4.75) 
Convergence (|4.74p is proved in Step 1 and (|4.75p follows as in the proof of Theorem [ 



Step 4 (Convergence of the load term). The statement follows directly from the condition (|4.40p . 
Definition 13.81 and Proposition I A. 41 

Step 5 (Convergence of infima). For brevity we set 

:= liminf inf A^((i3) — / w- Udx] . 
eo:= min ( A'"s<=(50, gi) _ £0( o i) ) 

We have to prove that lim = eg. Note that the inequality limsup^ < eg follows directly from the 

the existence of recovery sequences, and the property that the limiting energy admits a minimiser. 
Thus, it suffices to show that 

liminfej>eo. (4.76) 
Let ((^e) C iJp o(^) ^ sequence of (almost) minimisers, i.e. 



lim sup 



(4.77) 



and {g^,gl) *^ (pe the associated splitting. By Theorem 13.41 and Proposition 13.101 we can pass to 

a subsequence (not relabeled) such that (p^ (.9"i3^) for some {g'^,g^) G Qr.oi^l'xY). Using the 
decomposition estimate of Lemma 14.101 the identity from Lemma 14.111 and by appealing to the 
definition of and we rewrite the total energy as 

UVe) - f ^e- le dx = Ag (uM)) + Kigl) - ff ffs ' dy dx - f gl ■ ll dx + 0(£^). 

Setting := Ueide) can express the right-hand side of the last identity by means of the limiting 
energy: 



A^icp,)- I ip,-kdx = A'-s°(G„5^)-/:"(G„gi) 
+ {Aligl)~Alig')) 



Ge-{r -r;)dxdy (4.78) 

QxYo 



{g^ -l^ -gl-lDdx + Oie-^). (4.79) 



Note that the "error terms" (|4.78p and (|4.79p vanish as e — !• 0, as can be seen from ProDOsition l3.101 bv 
taking into account the assumption (|4.401) . Since g]. g^ in H^{VL), we have Mvami A\{g].) > Al{g^) 
and hence arrive at the lim inf estimate 



lim inf 



nf (^A,{^e) ~ l/^- dx^ > liminf (A'-s^(Ge, g^) - C°{G„ g^)) > cq, 

where in order to get the last inequality we employed the minimality property of eg. Since the 
right-hand side is independent of the specific choice of the subsequence, (|4.76p now follows from 

627D. □ 
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5 A version of the geometric rigidity estimate for perforated 
domains 



A central ingredient in the proof of the a priori estimate of Theorem 13.41 is a geometric rigidity 
estimate for the perforated domain ^\ , which allows one to control the (fi^ ) distance to a single 
rotation of Vu, u £ H^{fl), by the integral J^^^ dist^(Vii, SO((i)) da; with a constant independent 

of £. The geometric rigidity estimate for bounded Lipschitz domains (see Theorem 15. ip has been 
developed in a seminal work by Friesecke, James and Miillcr in |23j . where the following statement 
is proved. 

Theorem 5.1 (Geometric rigidity estimate, see [23] ). Let be an open, bounded Lipschitz domain 
in R'', d > 2. There exists a constant C{fl) with the following property: for each v G H^{il) there is 
a rotation R G SO{d) such that 

[ \Vv{x)- R\^dx <C{n) [ dist2(Vu(x),S0(d))dx. 
Jn Jn 

Moreover, the constant C(f2) is invariant under uniform scaling ofQ. 

In [23j the above estimate leads to a rigorous derivation of the nonlinear bending theory for 
elastic plates from three-dimensional finite elasticity. Apart from this, it is also the key ingredient 
in the rigorous justification of linearised elasticity from finite elasticity as established in [20]. In 
this section we present a simple argument that shows that the rigidity constant associated to the 
perforated domain 0,1 can be chosen independent of e. We would like to remark that, independently 
of our work, a related result for domains with a finite number of circular e-holcs has been recently 
established in an unpublished work by Scardia and Zcppieri |35| in their analysis of dislocation 
energies. 

Proposition 5.2. Let 0, Yi and fll be as assumed in Sectionl^ There exists a positive constant 
C(f2, Yi) that depends on the geometry of O and Yi only, such that the following holds: for all e > 
and V G H^(Q,l) there exist a harmonic extension v G H^{Q,) of v (i.e. u = w m fif and v satisfies 
(|3.9p ) and a rotation R G SO{d) such that 

||dist(ViJ,SO(d))||^.(,,) < C(l],yi)||dist(Vi>,SO(d))||^.(^,), (5.80) 

\\Vv-R\\mni) < Cin,Yi)\\dist{Vv,SO{d))\\mni)- (5.81) 

Proof. Since (|5.8ip is a direct consequence of (|5.80p and Theorem 15. 1[ it suffices to prove (|5.80p . 
To this end let v denote the harmonic extension of v to 17 as introduced in Lemma 13.21 Define the 
following sets 

Z,:={^eeZ'^:^ + eYiCnl}, 17^ := |J (^ + eF). (5.82) 

By the definition of fll, the functions v and v are identical on \ fJ^. Thus, it suffices to argue that 

/ dist^(Vw(x),SO(d))da; < C / dist^{S/v{x),SO{d)) dx. (5.83) 
Jn^ Jn^nni 

For the argument, we introduce the re-scaled functions U^{y) := e^^v{£^ + ey) defined for y G y and 
^ G Z^. By construction, we have WU^{y) = Wv{(, + ey); therefore, estimate (|5.83p follows, provided 

[ dist^{VU^{y),SO{d))dx<C [ dist^{VU^{y),SO{d)) dx (5.84) 

Jy JYi 

holds for all ^ G and a constant only depending on Yi. Indeed, this is the case: let R^ be a 
minimiser of jy^ l^^c(y) ~ -^C P SO((i). Since y i— > U^{y) — R^y is harmonic in Fqj we find that 

/ \VU^-R^\^dy<C [ \VU^~R^\^dy. 
Jy Jyi 
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By the geometric rigidity estimate, Proposition 15. 2[ the right-hand side of the previous hne is 
estimated by the right-side of ()5.84|) (up to a constant that only depends on Yi). Moreover, the 
left-hand side of the previous inequality is clearly an upper bound for the left-hand side of (|5.84p . 
In conclusion, (|5.84p follows. □ 

In situations when the boundary data of the deformation v in Proposition 15.21 are known, one 
can gain further information on the rotation R. Suppose for a moment that v{x) := x + f{x) with 
(fi G Hq{H). Then by partial integration we have 

/ \\7v{x) ~ R\^dx = I \^ip{x) - [R- dx ^ I |V^(a;)|2dx+|r2||i?-/p, 
Jn Jn Jn 

which, together with geometric rigidity, implies 

/ \\/v{x) ~ I\^ dx < ( \\/v{x)-R\^dx<C(Ji) [ dist2(Vi;,S0(d))da;. 
Jn Jn Jn 

In the situation above we could eliminate arbitrary rotations, by appealing to the Dirichlet bound- 
ary condition. A refined version of this argument was given in ^0] , where also non- homogeneous 
boundary data arc considered. For us it is sufficient to borrow a technical statement from [20j : 

Lemma 5.3 (c/. Lemma 3.3 in [20]). Let il be an open and bounded Lipschitz domain in Mf^ and F a 
measurable subset of dil with positive {n — 1)- dimensional Hausdorff measure "H""^. For F e R'^'^'' 
we define 

\F\l := min / iFx - 5P dH"" Vx). 

bGK'' Jy 

Then there exists a positive constant C such that 

\F\^ < C|F|p 

for all matrices F that belong to I ~ SO(ri). 

Now, we can easily deduce the following result: 

Proposition 5.4. Let $7 and F as assumed in Section\^ Then there exists a constant C = C(r2,F) 
such that for all u S iJ^(f2) and R € SO(ri) with u(x)\r — Rx the estimate 

II - R\\mn) < C ||dist(Vw, S0(7i))||^.(f,) 

holds. 



Remark 5.5. Note, that the statement of Proposition \5.4\ is valid also for fl^ with a constant 
independent of e. It is easy to see by applying consequently Provosition \5.2\ and \5.4\ We make use 
of this observation in the argument of Section\^ above. 



Proof of Proposition \5 .4\ According to Theorem 15.11 there exists a rotation P G S0(7i) with 

II - P|li.(o) < cm ||dist(Vu, SO(n))||i.(^) . 

Thus, 

II - Rfmn) < ^C{n) \\dist{Vu, SO(n))||i.(j,) + 2 \n\ \P - R\' . (5.85) 

Define the map v{x) := u{x) — Px + Vfi where vn € R'' is chosen in such a way that v has zero 
mean value over n. Then by the Poincare inequality we have 

\Mm{n) < c||Vu||^2(n) < cC(0) ||dist(Vu, SO(7i))||^2(a) , 

since Vw = Vu — P. Within this proof, c denotes a constant that may vary from line to line, but 
can be chosen in such a way that it depends on fl and F only. 
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On the other hand, we have v{x) = u{x) — Px + vq = {R — P)x + wji on F. Recall the definition 
of the semi-norm |-|p from Lemma l5.3l By construction we see that 

|i?-P|r< / \{R-P)x + vn\'dn-\x)< [ \v{x)\'dn-\x)<c\\v\\l,^^y 
Jr Jon 

In the last estimate, the right-hand side is controlled by ||dist(VM, S0(n))||^2(Q) and by virtue of 
Lemma [5. 3 [ the left-hand side controls the distance \R — P\^ . This implies that 

|i?-P|2<c||dist(Vii,S0(n))||i.(^) 

and, combining this with ()5.85p . the statement of the theorem now follows immediately. □ 

A Appendix: Two-scale convergence 

In this section we gather basic definitions and properties of two-scale convergence (resp. periodic 
unfolding) following the lines in [52]. Let J7 C M'' be open, bounded and satisfy \dn\ = 0. Following 
[Tl [32] a bounded sequence (ue) C L^(ri) is two-scale convergent to a function u £ L'^{ilxY), if for 
all test functions i/' € L^{^; C^{Y)) we have 

lim / Ui,{x)'ip{x,x/e)dx ^ // u{x,y)ip{x,y)dxdy. 
^^°Jn JJnxY 

In the present paper it is more convenient to work with the equivalent definition of two-scale con- 
vergence via periodic unfolding (see [T^], [37], [IH])- On M.'^ we define the mappings [•] and {•} such 
that 

[■] : ^ Z'' , {•}:M'^^y, x = [x] + {x} for ah a; e M''. 

Definition A.l (Unfolding operator). Let il C M.'^ be open and u : ^ M. a Lehesgue-measurahle 
function. The periodic unfolding of u (on scale e > is defined via 

Ue{u) : E^xF^M, {x,y) ^ u* {e[x/ e] e{x/ e}) 

where u* denotes the extension of u to by zero. For vector valued maps the periodic unfolding is 
defined componentwise. 

The following statements are well-known and can be found in [16], [39] . 

Proposition A. 2. 1. (Isometry property). For all s > the mapping if i— )■ ihti^ip) defines a 
(non-surjective) linear isometry from L?'{Vl) to L?'{W^'xY). 

2. (Product rule). Ifu,v E L'^{n), thenU,{uv) ^Ue{u)Ue{v) e L'^{R'^xY). 

3. (Integral identity). Let f : YxM.^ — > [0,oo] he a Borel map. Then for all measurable maps 
if : Vl ^M!' one has 

f{{xle},ip{x))dx ^ Ij f{y,Uei(p){x,y))dydx. 

J JR''xY 

Definition A. 3. Let (u^) C L^{n), u G L'^{VlxY), and let u* denote the extension of u to W^xY 
by zero. 

1. We say that (ue) strongly two-scale converges to u and write "u^ ^ u in L'^iyixY)", if 
Ueiue) u* strongly in L^{R'''xY). 

2. We say that (u^) weakly two-scale converges to u & L'^{QxY) and write "u^ u in L^{VLxY) " , 
ifU^iue) u* weakly m L'^{R'^xY). 
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Proposition A. 4. 1. If ^ u in L'^{ftxY), then \ \ue\\L^^n) is bounded for e — ;> 0. 

2. If Ug u in L'^{D.xY) then for any Y -periodic a £ L°° (R) one has a[x/e)ue{x) a[y)u[x,y) 
in L^{n X Y). 

3. We have ^ u in L'^{VlxY), if and only if u^ ^ u in L'^{{lxY) and ||u£||L2(n) ~^ 

\\u\\L^(ny.Y)- 

4- // limsupg_).Q ||we||i2(-Q-) < oo, then, up to a subsequence, ^ u in L'^{VlxY) for a map 
u £ L^{nxY). 

5. If Ug u, Ve ^ V, in L'^(n X Y), then J^UeV^dx^ /Q^yWwdxdy. 
Proposition A. 5. Let (u^) be a sequence in II^{il) and u £ 

1. If Ue ^ u in H^{n), then, up to a subsequence, Vug Vu + ^yip in L'^{^xY) for a map 
^£L\n-Hl{Y)). 

2. If We and eVii^ are bounded in L^(f2) as e — > 0, then, up to a subsequence, Uq in 
L^iflxY) and Vu^ ^ VyUQ in L^{D,xY) for a map uq £ L'^ (fl; H^{Y)) . 

3. If Ue £ HQ{rt'^) for all e > and limsupg_^Q ||Vue||^2(f20') < oo, then, up to a subsequence, 
ll"e||L2(n) = 0{e) and Vu. ^ Vyip{x,y) for a map -0 £ L^{fl]H^{Ya)). 

4. If Us £ i?o (ri") for all e > and lim supg_j.Q I jeVuj I |^2^f^o) < 00, then, up to a subsequence, 
Ue Uq in L'^l^lxY) and Vu^ ^ V^itQ in L'^l^lxY) for a map uq £ L'^ (fl; H^{Y)) . 

Lemma A. 6 (Attouch's diagonalisation, see [2]). Assume e,6 > 0, h : [0, 00)^ — > [— 00,+cxj]. Then 
there is a mapping e 1— > (5(e) such that 

lim (5(e) —>■ and lim sup /i(e, (5(e)) < lim sup lim sup /i(e, (5). 
In particular, if h > and limA-_i.o liine->o ^{e, S) = 0, then lim£_i.o ^{e, (5(e)) ~ 0. 
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